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ABSTRACT 


The Modular Modeling System (MMS) was developed by the 
Electric Power Research Institute and the Babcock and Wilcox 
Company for the study of the thermal-hydraulic performance 
of pressurized water reactor plants at steady state and 
during some slowly varying transients. In order to 
determine the limits of transient severity that this code 
can endure, two experiments conducted at the Loss-of-Fluid 
Test (LOFT) facility at the National Reactor Testing Station 
in Idaho were simulated using the IBM/CMS computer of The 
Pennsylvania State University. 

These experiments are a small break loss of coolant 
accident (Experiment L3-5) and an excessive steam generator 
load increase (Experiment L6-3). In the case of the former, 
the Modular Modeling System failed to accurately predict the 
performance of the LOFT facility. The MMS was, however, 
successful in predicting the significant thermal-hydraulic 
parameters of Experiment L6-3. The MMS predictions of the 
LOFT facility's performance during this experiment were more 
accurate than those of the more sophisticated RETRAN code. 

This success validates the MMS' ability to predict the 
performance of nuclear power plants that are scaled much 


smaller than central station plants. 
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Chapter l 


INTRODUCTION 


Pressurized water reactor safety has been a key concern 
of environmentalists, various governmental agencies, 
electrical utilities, and reactor vendors since the 
beginning of the nuclear energy era in the 1950's. Today, 
with major nuclear accidents having occurred in both the 
United States and the Soviet Union, reactor safety is again 
at the forefront of issues driving the nuclear industry. 

Fundamental to the issue of reactor safety is the ability 
to predict the performance of key plant operating parameters 
during both slowly varying transients and catastrophic 
accidents. These parameters can be separated into two 
general areas: core neutronics and system wide thermal- 
hydraulic characteristics. Thermal-hydraulics, in 
particular, is the focus of interest when a loss of coolant 
occurs from the primary volume of a pressurized water 
reactor (PWR). Thermal-hydraulics is the study of the the 
heat transfer and heat transport properties of fluid 
systems, and has been the focus of study of physicists and 
engineers for many years. These professionals have 
successfully reduced both core neutronics and plant thermal- 
hydraulics to a series of fundamental equations. Accident 


prediction work is based on these equations. 





Although these fundamental equations are relatively few 
and simple, their use in reactor plant safety studies 
requires complex numerical methods for solution of the many 
resulting simultaneous equations. These complex numerical 
methods, in turn, require the use of computer codes written 
in today's modern languages to effect predictions within a 
reasonable amount of time, and with reasonable accuracy. 
The Babcock and Wilcox Company, in conjunction with the 
Electric Power Research Institute, has developed such a 
computer code, the Modular Modeling System (MMS). This code 
was not designed specifically to predict the performance of 
pressurized water reactors under accident conditions. 
Rather, it is intended for the early design stages and gross 
predictions of any large electrical generating plant, be it 
fired by conventional fossil fuels or a nuclear reactor. 
However, this code does have some features that make it very 
desirable for reactor plant accident studies. These 
features are its modularity concept and its fast execution 
time. 

There are other computer codes that are suitable for 
reactor plant accident studies. Those that were designed 
specifically for reactor plant use include the many versions 
of RETRAN (developed by the Electrical Power Research 
Institute) and RELAP (Idaho National Engineering 
Laboratory), TRAC, and SIMMER (Los Alamos National 


Laboratory). For a variety of reasons, these codes are more 





Suitable for the in-depth investigation of reactor plant 
performance than is the Modular Modeling System. The 
expense required by their use, however, demands a cheaper 
alternative. Hence the Modular Modeling System became a 
candidate for the prediction of gross reactor system 
performance. 

Because the MMS was designed for examining power plant 
performance in the steady states in addition to some small 
transients, it is logical to attempt to extend its use to 
include more severe reactor plant accidents that do not 
rapidly change the physical states being modeled by the 
Simulation language. 

In order to determine if the MMS indeed can be relied 
upon for accident analysis, its performance in predicting 
key reactor plant parameters in an actual accident needs to 
be known. Many accidents and unuSual occurrences have 
happpened during normal operation of reactor plants in the 
United States. However, the instrumentation systems of the 
typical utility's power station are not designed to record 
the immense amount of data that is required for in-depth 
accident studies. 

Other facilities have been established for the specific 
purpose of accident analysis, including the initiation of 
accidents on actual reactor plants. One such reactor plant 
is the Loss-of-Fluid Test (LOFT) facility at the National 


Reactor Testing Station near Idaho Falls, Idaho. This fully 





Operational reactor plant was designed and built 
specifically to provide a test bed upon which actual reactor 
Plant accidents could be initiated and studied. These 
accidents range from small break loss of coolant transients 
without the nuclear core installed, to large break events 
that began with the core at 100 power. Because of the 
deliberate nature of the tests at the LOFT facility, 
instrumentation was in place that provided accurate records 
of the thermal-hydraulic characteristics of the plant 
throughout the various accidents. These records are the 
basis upon which the MMS can be evaluated as a useful code 
in accident studies. 

Because the Modular Modeling System was never intended to 
model such severe transients as would occur ina large break 
loss of coolant scenario, an attempt to determine its 
ability to be used in large-scale accident studies would be 
doomed to failure. Instead, a smaller scale accident in 
which the parameter changes occur more slowly is the choice 
upon which such a determination could be made. Since many 
such accidents were performed at the LOFT, a selection of 
two of these experiments was used as the basis for this 


validation study. 





Chapter 2 


THE MODULAR MODELING SYSTEM 


2 1 Objective 


The Modular Modeling System was developed to provide 
an easy-to-use, flexible, economical, and accurate 
systems analysis code that can be used for 
Simulating and analyzing the dynamic performance 
of nuclear .. . power plants. 
This effort was conducted primarily by the Electric Power 
Research Institute of Palo Alto, California. Other 
participants include the Babcock and Wilcox Company of 
Lynchburg, Virginia, and the Bechtel Group of San Francisco. 
Specifically, this system is intended for 
-specification, selection and integration of plant 
components 
-design and checkout of control systems 
-rapid simulation to expedite plant commissioning 
-best estimate plant sgfety analysis 
-procedure evaluation. 
These objectives were designed to fit the constraints of 
minimized computation time, minimized time for model 
development, and a reasonable amount of confidence in the 


generated results. The primary characteristic of the MMS 


included to meet these objectives is the modularity concept. 


+ R. R. Dixon, S. W. W. Shor, and Lance P. Smith, The 


Modular Modeling System (MMS): A Code for the Dynamic 
Simulation of Fossil and Nuclear Power Plants (Palo Alto, 


California: Electric Power Research Institute, 1983), I, v. 


- Dixon, Shor, and Smith, I, 1-l. 





2-2 Modularity Concept 


The current library of the MMS includes 64 modules. 
Most are designed to represent a typical component of an 
electric power generating station. These components range 
from the hydro-mechanical, represented by pipes and valves, 
to the electro-mechanical, represented by the on-off 
controller. The modules are divided into six basic groups 
which depend on power plant type. The fossil group includes 
component modules such as oil and coal fired boilers. The 
nuclear group includes reactors and steam generators. The 
controller group has_ proportional-integral signal 
generators, and the extended range group includes newer 
modules that allow for two-phase fluid flow. The balance- 
of-plant-component group includes pumps, pipes, and valves; 
components found at all power stations. Finally, the 
general group includes connections and junctions. This 
group iS unique in that its modules do not always represent 
actual physical plant items, but are necessary to meet the 
connection requirements explained in the next sections. 

The MMS modules have been designed to be joined together 
in a configuration which the user determines to best 
represent the actual physical system he/she wishes to 
Simulate. In this joint configuration, all the physical 
properties the user wishes to calculate would then be 


determined by a FORTRAN computer program. This program, 





generated upon completion of the module joining process, 
will perform its calculations based on the physical 


properties discussed in Chapter l. 


gree. Module Joining 


Each of the modules which simulates the containment of a 
Flowing or static fluid is called a hydro-mechanical module. 
These modules are further sub-classified as resistive, 
storage, resistive-storage, or storage-resistive components. 
The fluid in these components 1s usually, but not 
necessarily, restricted to a single physical phase. 
Usually, either a liquid or a vapor phase is modeled. The 
vapor can be either saturated or superheated. 

These modules must be joined in such a manner that the 
fluid "flows" alternately between resistive and storage 
nodes. The modules were designed so that adjacent modules 
do not solve for the same physical property, for example 
pressure, and no property is left undefined in a closed loop 
Flow system. In system models that have open ended flows, 
the properties at the flow boundaries are maintained by user 


selected boundary conditions. 


2 Ja ae Resistive Modules 


PeReodule which uses a pressure drop as the basis of its 


Flow calculations is described as purely resistive. A 





typical example is a simple pipe. The arrow convention in 
Figure 1 shows that the inlet and outlet flow stream 
pressures, P, must be supplied as numerical Ayanayute EG 
resistive modules. The inlet flow properties are shown at 
the left of the figure, and the outlet flow properties at 


Ewesraght. 


NZSISA Ve 


P 
W: MODULE 
h (PIPER) 





Figure 1: Resistive Module 


The other input required is the inlet enthalpy, h. The 
arrowS pointing out of the module indicate the calculated 
values delivered to the adjacent modules. These are the 
mass flow rates and the outlet flow enthalpy. Since 
resistive modules have no storage volume, the mass flow rate 


into the module must equal that flowing out. 





2) Storage Modules 


The continulty equation, 
ole [We - Ww Eqn. 2-1 


is used as the basis of the storage modules. As Equation 
2-1 and Figure 2 show, the inlet and outlet mass flow rates 


are needed as module input. 





Figure 2: Storage Module 


The primary output produced are, then, the inlet and outlet 
pressures. Of course, the other variables in the equation 
are assumed known. Such parameters as a tank volume, V, 

would be provided as input to the program. CONNI is the MMS 


name for a connective module. 


Bae. l.3 Storage-Resistive and Resistive-Storage Modules 


Modules that determine both pressure drops and flow rates 
combine the equations used in the purely resistive and 


storage modules into a single module. Such modules are 





10 


designated as resistive-storage or storage-resistive. A 
resistive-storage component is shown in Figure 3. Note that 
the inlet flow pressure, at the left of the figure, is by 
arrow convention a module input. This indicates the flow 
encounters the "resistive" section of the module first, and 
then the storage section. Simply reversing the P and W 
arrows would make this component into a storage-resistive 


type. UTSGR is a U-tube steam generator module. 


Pei aiye = 
STORAGE 
MODULE 
(UTSGR) 






= 16 


>- 


Figure 3: Resistive-Storage Module 


eid. 2 Connectivity 


Meeting the requirement that modules be joined in 
alternating resistive/storage fashion seems an easy thing to 
do. However, computing costs are very sensitive to the 
modules selected for use in the MMS, and hence the user must 
be very careful to select those modules which suit his/her 
individual requirements. Because the optimum configuration 


may require two resistive modules adjacent to one another, 





agi 


seemingly violating the alternate resistive/storage 
requirement, still other modules were developed. These are 
appropriately termed connective modules. An example is a 
10" pipe which is connected and constricted down to a 7" 
pipe. The connective module simply determines the missing 
pressure between the two resistive pipe modules. Since 
there is no envisioned configuration where two storage 
modules would be adjacent, there is no "connective" module 
of the resistive type to join them. All tanks, in reality, 
must be joined by a pipe (resistive module), no matter how 
short. 

It is not required that the modules be connected in 
Closed loop fashion. An open loop, or a closed loop with 
some inlet and outlet connections, iS permissible. The 
models developed for this study are of the latter type. Any 
connections left open must be modeled by boundary 


Gonaitions. 


Ze 3 Physical Models 


The physical phenomena upon which the MMS is based are 
listed in Table 1. These phenomena include some of the 
basic physical laws of thermal-hydraulics and heat transfer. 
All are applied to the appropriate modules, treating each 


module aS a separate control volume. 





Table 1: Basis of the Modular Modeling System 


Phenomena Basic Governing Equation 
Conservation of Mass Jes stovi) 
dt Xj 


Conservation of Energy d{oe} _ -dfoevi) — a{pvy] 
at 


0X; 9X; 


+q -W- 9(dccV; 


OX; 

Conservation of 3 Vie 
; OVale so. | OMI dp 
Linear Momentum afew) Pear id ae 
! 1 


— 9-9 0j; - og{sin 8) 


OX; 
Radiant Heat Transfer q = UA(AT) 
Convective Heat Transfer Dittus-Boelter Equation 


(other heat transfer 
equations are listed 
in reference [5]) 


Viscous Shear Losses Ao meen Vi 


KZ 





aS. 


The conservation laws are integrated around the surface 
of each module's control volume as the user directs. For 
example, he/she can include or ignore heat losses from a 
nuclear station's primary plant piping. To ignore these 
losses the MMS would simply integrate using the fluid 
boundaries into and out of a pipe module. To include these 
losses the pipe walls themselves would become part of the 
integration boundaries when performing the energy 


conservation calculations. 


ees. 1 Liebnitz's Rule 


In order to apply these physical laws to the thermal- 
hydraulic characteristics of an actual power generating 
system, a set of ordinary differential equations must be 
developed by the modules. This process makes use of 
Liebnitz's Rule to arrive at the integral form of the basic 


equation. 


fv 30 14 Jloav—f. ovg.dA Eqn. 2-2 
dt dt 


The left side of the Liebnitz Rule equation is the integral, 
over a given control volume, of the rate of change of a 
certain thermodynamic property. The right side of this 
equation is the form which must be set equal to the three 
phenomena of the conservation laws. 

To illustrate the application of this rule, it will be 


applied to one of the conservation laws. The simplest 
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example is that of the conservation of mass. Converting 


Peuation 2-1 to integral form yields 


iv go dV _ —\y a(ovi) dv 


Bane 2-35 
at aX; 
which using the Liebnitz notation becomes 
d WE =- Vi V 
ee. end JsoVs dA J d\0 d Eqn. 2-4 


The right or divergence term is converted to a surface 


integral using the divergence theorem, 





1 ON me fo Vid A Eqn. 2=5 


thus changing Equation 2-4 to 


4 Jyodv= fs ovg dk - fs waa Eqn. 2-6 
dt 


The left side of this equation is the change in the 
instantaneous mass in the control volume, dM/dt, and the 


surface integrals on the right reduce to 


RHS = AgdaVSe * AIDIVS * AePeVe - AIIM , ei I 
Hence the conservation of mass can be written as 
4 / 
dM = AgdalVE - Vge) - Arpi(V) ~ Vsi) Eqn. 2-8 
dt 
which contains no partial derivatives. 
Similar developments, described in reference [5], are 
applied to the Conservation of Energy and the Conservation 


of Linear Momentum. These laws describe the "“bulk" 
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properties within a control volume, properties that are 
termed extensive. Their equations alone are not enough to 


provide a solvable set of differential equations. 


wes! Intensive Properties 


The remaining properties of the fluid needed to form a 
complete set of differential equations are labeled 
"intensive" properties. These include the specific density, 
temperature, and internal energy, and may have different 


values at the inlet and outlet boundaries. 


ma. 2.1 Extensive/Intensive Property Relationship 


To show the relationship between the intensive and 
extensive properties of each module, the average fluid 
enthalpy (intensive) and total internal energy (extensive) 
can be examined. In a control volume which has no heat 
transfer across its boundaries, the change in total internal 
energy is 

dU = Wehe ~ wyhy Eqn. 2-9 

dt 
when there 1s no work done by friction. Using h as the 
average enthalpy within the the control volume, and M as the 


total fluid mass, 


dh = Wane > why - (dM/dt)n 7 
i M imo fole A IE0) 
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is the change in the specific enthalpy per unit time. 
Assuming that the change in h as the fluid flows through the 
module is linear, and that the inlet enthalpy is provided by 
the upstream module or is a boundary condition, then the 
outlet enthalpy 1s easily computed. Any instantaneous 
changes in the inlet property are mitigated at the outlet by 
assuming the derivative of the property leaving. the node is 
equal to the derivative of the average value. 

Forming a complete set of differential equations for each 
module is, then, a two step process. First, the extensive 
properties are used to determine average values of the 
intensive properties, and then these average values are used 
in the determination of the outlet flow intensive 
properties. These outlet values are then used as the input 
to the adjacent modules. For example, the derivative of the 
Specific energy leaving the control volume becomes 

dU) _ 1 { Wehe - wih) tq -W - Eom. 2-11 


dt ov a . 
p(dV/dt) - Up(dV/dt) - UV(de;/at)} 


after conversion to the intensive form. 


25 Sa The Complete Set of Thermal-Hydraulic Equations 


The MMS uses enthalpy and pressure as the system states 
upon which its solutions to the sets of differential 
equations are based. The developers selected these 


properties because they are commonly used in many 
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engineering applications and to minimize numerical stability 
problems known to occur when basing models on density. 
Selection of the equations used to solve for 
thermodynamic properties depends on one final factor, the 
number of fluid phases present within the control volume. 
In the case of a Single phase, the equations are those 
developed above, one each for the conservation of mass and 
the conservation of energy. There are two differential 


equations and two unknowns, enthalpy and pressure: 


dp} — 1 {1 (We - wy - p(dV/dt)] - a_(dh/at) } Eqn. 2-12 
dt a V 

p 
dhy — _1 (Wehe - why + q - Ws - Eqn. 2-13 
dt oV 


oh(dV/dt) - hV(dp;/dt) + V(dp;/dt)) , 


In these equations, the properties other than enthalpy and 
pressure are determined by FORTRAN steam property 
Subroutines. The heat added to a control volume, gq, 
requires a separate calculation using the heat transfer 
equations of Table 1. This separate calculation does not 
have any unknown values, but rather depends upon results 
from the previous time step, or boundary values. 

Two-phase control volumes require the use of the "drift- 
Elux’ concept. "Using the drift-flux model and the 
assumption of inter-facial equilibrium between the steam and 


liquid phases, the separate phase equations can be combined 
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into a single set."2 The use of two phases in a single 
volume introduces such variables as the quality, void 
Fraction, and drift velocity. As the term "Separate phase 
equations" implies, the approach is to treat each phase as a 
separate control volume within each module, and introduce 
other unknown terms that account for the "drifting" of mass 
and energy across the boundary from one phase to the other. 


The equations used are 


ap; dh; + d0; dp; l 
efile EE = es Ean. 2-14 
dh; dt pj dt es v2) 2 : 


for the conservation of mass and 


(pj + Hj(Apj/dhj))hj + ((9p;/ 9p;)h; -— 1/Je)p; = 
EGinleeea— lo 
{f7w 1H - fow2H2 + (1 - fy)wyH2 - (a - fo)w2HRe2 + a}/V 


+ ACf1V 421 - foVoz2 + (1 - £1)V122 - (1 - f2)V2zR}/V 


for the conservation of energy. 
In a two-phase storage module, the level of the liquid is 
determined using the average void fraction and the densities 


of the liquid and vapor. 


mes. 3 Reactor Kinetics 


The reactor modules available in the MMS use the point 
kinetics equations to calculate reactor power. These 


equations are 


3 Dixon, Shor, and Smith, I, 3-16. 





ro 


oe (1 = Bel; +) A;Cjj + Oe Meare) = May | EG ie 2 1.6 
dt Af iy 
and 
<i Bj (Ni/A) = AGCj. FE Gileceeee ee 7 
t 


The number of groups of delayed neutron precursors varies 
with the module used. The third term on the RHS of Equation 
2-16 accounts for the migration of neutrons from one node to 
another in the multimode modules. 

Once the number of fissions occurring in an integration 
time interval has been calculated, the heat produced in the 
fuel is determined. This heat is then the basis for the 
change in enthalpy in each node uSing the system of 


differential equations described earlier in this chapter. 


ees. 4 Other Physical Processes 


The MMS accounts for the processes of viscous shear 
losses and heat transfer according to the basic equations 
listed in Table 1. The process of transport delay is 
accounted for by simple memory delays or by the use of 


multiple nodes (modules) to represent piping runs. 
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a il Integration 


The MMS is based squarely upon the use of numerical 
integration techniques to solve the differential equations 
explained above. Each model creates a matrix of equations, 
the size of which is roughly proportional to the number of 
modules in the model. This matrix is solved for the unknown 
enthalpies and pressures at a fixed or variable time step. 
Once these new values have been obtained, calculations using 
regular FORTRAN statements and steam property subroutines 
are made to determine the new values of any other variables 
that may have changed. Such variables include tank levels, 
temperatures, and the extensive thermodynamic properties of 
each module. 

Many algorithms have been developed through the years 
that can be applied to solving these sets of equations. The 
MMS, because it primarily calculates fluid mechanics and 
energy equations, generates a certain range of time 
constants. The inverse of these time constants are the 
System eigenvalues. It is the range of these eigenvalues 
that determines the optimum algorithm used to solve the 
differential equations. The algorithm used in this study is 


the Gear's Stiff. 
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The Gear's Stiff algorithm is a "variable step, variable 
order integration routine that is self-initializing."4 This 
algorithm attempts to keep the size of each derivative in 
each time step below a preset value. If the size exceeds 
the set value, the time step is reduced and the set of 
equations is solved again. It also determines if the 
derivatives are approaching zero, and will increase the size 
of the time step if possible to minimize computation costs. 

The time constants of power plants are in the 0.1 to 100 
second range. The constants associated with the continuity 
equations are on the order of 0.01 seconds. This shorter 
time constant "stiffens" the process of solving the 
differential equations by implying the use of steps shorter 
than the user desires. Hence the term "stiff" is applied to 
the overall system. Simple Euler type algorithms must use 
the shortest time constant present as the time interval of 
the system solution and so would require about 10 iterations 
to advance even the smallest system time interval. Stiff 
algorithms make the assumption that the system's largest 
eigenvalues (shortest time constants) are always stable 
regardless of step size. This assumption can be made 
because the continuity equations are at quasi-steady state 


compared to the system time constants. 


4 Mitchell and Gauthier, Assoc., Inc., Advanced Continuous 


Simulation Language (ACSL) User Guide/Reference Manual (New 
MoEke Mitchell amd Gauthier, Assoc., Inc., 1981), p. 4-3. 
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In its current version, the MMS offers two "stiff" 
algorithms. The Gear's Stiff was selected because it will 
change the time step used in the solution many orders of 


magnitude, and so may possibly minimize CPU time. 


hs) Control Modules 


Modules which perform control functions use simple 
comparative calculations to determine their output signals. 
For example, the on-off controller is in reality a simple 
Switch with variable on and off setpoints. Other more 
complex controllers can have variable inputs and outputs 
which are determined by a series of polynomial equations. 
In this study, the most complex controllers used were the 
on-off type. The more complex versions were not required. 
Control modules do not contribute differential equations to 


the overall model's equation matrix. 


2.6 Model Generation and Execution 


To form and execute a complete MMS model, the designer 
follows the steps of Table 2. 
Figure 4 shows the flow process of steps 6 through 11 in 
Table 2, which are performed by the user's computer. At The 
Pennsylvania State University, the MMS and its host 
language, the Advanced Continuous Simulation Language 


(ACSL), are available only on the CMS system. 





1) 


2) 
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Table 2: Model Formulation Process 


Select the modules which best represent the 
system components needed to create a complete 
model. 

Collect the information needed to complete the 
applicable parameter sheets of reference [5]. 
Sources include plant data and operating logs, 
and vendor specifications and drawings. 
Determine the value of the initial operating 
parameters for use as boundary values. 

Write an ACSL program, using the required MMS 
syntax. 

Write an ACSL command file, using the necessary 
statements of reference [9]. 

Compile the ACSL program, locate and correct any 
errors discovered by the translator. 

Correct any errors discovered by the FORTRAN 
compiler. 

Create and load a FORTRAN module of the model. 
Execute the module through one iteration to 
determine if any time derivatives exist that 
exceed allowable error criteria, and make needed 
corrections. 

Execute the model to achieve steady state. 
Execute any transients of interest. 





Begure 4: 
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Chapter 3 


Abe LOSS OF -FLULD TEST SPACILITY 


Bei Objective 


The Loss-of-Fluid Test facility was designed and built to 
provide the United States with a capability to 
Simulate the major components and system responses 
of a commercial PWR during loss-of-coolant 
accidents (LOCAs) and during anticipated 
transients caused by abnormal PWR operations.? 
The facility is a fully operational 50 MW(t) pressurized 
water reactor plant designed to simulate the major primary 
System components of a commercial sized nuclear power 
generating station. This facility was erected at the Idaho 
National Engineering Laboratory near Idaho Falls, Idaho, in 
the early 1970's. 

In general, the intent in scaling the LOFT facility toa 
full sized PWR plant was to use the ratio of core power: 50 
MW(t) to 3000 MW (t). This ratio was "used as extensively 
as practical." 


The LOFT facility was subjected to many transients, 


ranging from small break "mini-blowdowns" without the 


2 Charles L. Nalezny, Summary of Nuclear Regulatory 
Commission's LOFT Program Experiments (Idaho Falls, Idaho: 
EG&G Idaho, Inc., 1983), p. 1-1. 


6 Douglas L. Reeder, LOFT System and Test Description 
(5.5-ft Nuclear Core 1 LOCEs) (Idaho Falls, Idaho: EG&G 


iano, Ine., 1978), p. 12. 
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reactor core installed to full-scale large break losses of 


coolant initiated at the reactor's maximum rated power. 


3.2 Primary Coolant System 


The primary coolant system removes the heat generated in 
the reactor core and transfers the heat to the steam 
generator, where it is passsed to the secondary coolant 
system. The primary coolant system also acts to contain any 
fission products that escape a fuel pin, and, with boron in 
solution, has a reactor control function. Nominal system 
pressure is 2250 psia, and rated flow at 100% power is 


3,780,000 lbm/hr. 


ee. L Intact Loop 


The LOFT facility has a two-loop primary system, shown in 
Figure 5. This figure shows that only one of the primary 
loops contains a steam generator and operating primary 
coolant pumps. This loop, called the intact loop, simulates 
three of the four loops of a Westinghouse Nuclear Steam 
Supply System (NSSS). It can be used to simulate any actual 
NSSS' coolant loops which have not been opened in a Loss of 
Coolant Accident (LOCA). 

Attached to the intact loop's hot leg (reactor outlet) is 
the pressurizer. This component's primary function is to 
maintain primary system pressure within the desired limits. 


The pressurizer is shown in Figure 6. 
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The steam generator is a vertical U-tube type (Figure 7), 
Similar to full-scale steam generators at actual power 
plants. 

The LOFT facility reactor coolant pumps are best 
described as the canned rotor, single stage centrifugal 
type, Similar to the pumps found at commercial nuclear power 
plants. A cutaway of one of these pumps is shown in Figure 


8. 


Be. 2 Broken Loop 


The facility's second loop is called the broken loop, and 
is used to simulate the large-scale losses of primary 
coolant that occur during a LOCA. This loop has no actual 
steam generator or coolant pumps. Instead, flow restricting 
devices called simulators create the pressure changes found 
across the steam generator and pumps in the intact loop. 
The other major components of this loop are the blowdown 
valves. These valves are of the quick-opening, 
hydraulically operated sleeve type. When opened, they act 
as the break location in a rapid loss of primary coolant 
experiment. The majority of the primary coolant system is 


constructed of 14" diameter stainless steel pipe. 
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633 Nuclear Reactor 


The nuclear reactor is the heat source of the LOFT 
facility NSSS. Like the reactors of central power stations, 
this reactor has a cylindrical single pass core, with the 
inlet and outlet pipe connections near the top of the 


pressure vessel. 


re S «1. Structural Components 


The reactor is contained within a 22 foot pressure vessel 
which is attached to the primary coolant loops as shown on 
the left side of Figure 9. This figure also shows the 
general internal arrangements of the reactor. The major 
components include the vessel's removable upper head and 
non-removable lower head, the core support barrel, flow 
Skirt, and the upper and lower core support structures. 
These support structures act to hold down the core against 
the pressures of the pasSSing primary coolant, maintain core 
and control rod alignment, and allow for thermal expansion 


and stresses. 


5) 6 4 Reactor Core 


The core contains the uranium fuel used to power the LOFT 
Facility's reactor. The core in place during the 


experiments of interest in this study was actually the 
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second used at the facility. This core is a modified 
cylindrical design with a length of 5.6 feet and a median 
diameter of 2.3 feet. 

The fuel assemblies in the core closely approximate the 
design used in actual power plants. There are five 15x15 
Pin square assemblies and four triangular corner assemblies 
that contain 12 pins on each side, making a total of 1300 
pins. The core is assembled as shown in Figure 10. The 
LOFT facility cores are rated at 2000 effective full power 
hours at 50 MW(t). 

Neutron generation rates in the core are controlled by 
four spider type control rods and the use of soluble boron 
in the primary coolant. The control rods are located in the 
four square fuel assemblies that surround the center 
assembly. Their neutron absorbing materials are silver, 
indium, and cadmium. Because the LOFT facility core is 
small compared with the commercial cores it simulates, there 
was no need to devise a bank rod control system. The 
control rods all move at the same speed and time during 


normal operations. 


3.4 Additional Primary Systems 


B44. 1 Emergency Core Coolant System (ECCS) 


In the transients examined in this study, changes in 


primary system pressure initiated emergency core coolant 
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system operation. This system was designed for plant 
protection by ensuring the core would remain covered with 
liquid coolant after any size piping break, and act as part 
of a long term shutdown cooling system. Performance of 
these functions was intended to be as close to the 
performance of an actual ECCS, with safety of the LOFT 


Facility being the overriding consideration. 


ee? .1.1 High Pressure Injection System (HPIS) 


The HPIS is designed to make up lost primary coolant 
during small and intermediate break events. This system has 
two positive displacement pumps and a nitrogen pressurized 
accumulator system to perform this function. Both act to 
inject borated makeup water into either the intact loop hot 


and cold legs, or the reactor vessel itself. 


ayo a Low Pressure Injection System (LPIS) 


The LPIS acts with the HPIS to mitigate the more severe 
losses of primary coolant. However, as its name implies, 
the primary system must be at a lower than normal pressure 
for the LPIS to operate. Such low pressures, if not 
accompanied by injection system operation, could lead to 
overheating of the reactor core in the worst case, or to the 
formation of non-condensible gases at the top of the primary 


coolant pump motor casings ina less severe case. 





3] 


The principal components of the LPIS are two single stage 
centrifugal pumps. These pumps take a suction on the same 
borated water storage tank as the HPIS pumps, but have a 
much higher capacity, 300 gpm vs. 17 gpm at nominal 
discharge pressure. The borated storage tank has a capacity 
of 24,000 gallons, roughly 12 times that of the primary 


coolant system. The ECCS is shown in Figure ll. 


3.4.2 Blowdown Suppression System 


The Blowdown Suppression System simulates’ the 
backpressure effects of the containment structure of an 
actual NSSS, and collects the discharges from the primary 
Piping during a fluid loss experiment. The major component 
is the blowdown suppression tank. This tank is a 
cylindrical vessel 38 feet long and 12 feet in diameter. It 
is connected to the broken loop by the blowdown suppression 
header and the quick-opening blowdown valves. There are 
other smaller connections to the primary piping, including 
One in particular to the intact loop cold leg that is 
important to this study. 

In order to cool the large amounts of very hot water 
discharged from the primary system during an accident 
Simulation, the blowdown suppression tank contains borated 
water at all times. The headers into the tank extend 
beneath the surface of this condensing pool. Additional 


cooling spray is also used during large break experiments to 
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condense the steam in the tank. This form of tank pressure 
control is used to best simulate a containment building. A 
diagram of the blowdown suppression system is shown in 


Figure 12. 


3.5 Secondary Coolant system 


ere 52.1 Steam Generator 


The heat delivered to the steam generator by the primary 
coolant system is transferred across the U-tube walls to the 
secondary side. This heat serves to change the feedwater 
from a subcooled liquid to a saturated liquid-vapor mixture. 
This mixture has a quality of about 25% when the system is 
Operated at rated power. The LOFT facility's steam 
generator uses two stages of driers to remove the moisture 
from the mixture: a Swirl separator at the top of the tube 
bundle shroud and a mist extractor just below the steam 
outlet. The steam moves vertically out the top of the 
generator's 23 foot length. The steam generator delivers a 
maximum of 220,500 lbm/hr of dry saturated steam at 808 psia 
to the condenser. This flow rate is controlled by the steam 
control valve, located between the steam generator and the 


condenser. 
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Saoe..2 Condenser 


The LOFT facility does not use its reactor generated 
power to produce any energy ina usable form. Instead, all 
the energy 1S removed from the steam by an air cooled 
condenser. This condenser consists of finned tubes across 
which forced air flows, moved by a set of variable pitch 
fans. There is only a small pressure drop across the fluid 
Side of the condenser, so the water arrives at the 
condensate receiver at very close to saturation temperature, 
520°F, unless the air flow is such that there is a larger 
than normal amount of subcooling. The cooling air is 


discharged directly to the atmosphere. 


oo. 3 Feedwater System 


From the condenser, the condensate flows into a 
cylindrical vessel called the receiver. This vessel acts as 
an expansion/contraction volume for the secondary coolant 
system. From here the condensate moves to a water cooled 
Ssubcooler. This subcooler is required to control the 
temperature and density of the water at the feedwater pump 
Suction, and so prevent feed pump cavitation. 

The feedwater pump is an electrically driven, multistage 
centrifugal pump which delivers the feedwater back to the 
steam generator via the feedwater regulating valve. fThe 


feedwater regulating valve is controlled by the steam 
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generator water level control system to match the steam flow 
rate with the feed flow rate and thus maintain the desired 
downcomer level in the steam generator. Also included as 
part of the feed system is an auxiliary feed pump which 
Supplies a small amount of makeup water to the steam 


generator when the main pump is inoperative. 


5. Sie 4 Design Theory 


It should be noted that the LOFT facility's secondary 
coolant system is markedly different from that of an actual 
nuclear generating station. The facility has no turbines, 
water cooled condensers, or electrical generators. fThis 
equipment is not needed because the primary intent of the 
facility is to examine the primary system performance under 
accident conditions. This unusual design makes modeling the 
facility's secondary system on a direct component by 
component basis impracticable because an air cooled 
condenser module has not been developed. However, some 
Substitutions and omissions can be made to create a suitable 
working model. Figure 13 is a simplified schematic of the 


secondary coolant system. 
8.6 Instrumentation 


Because the LOFT facility was designed for research, it 


was constructed with an extensive array of instrumentation 
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which uses fixed and movable detectors that monitor the 
parameters in Table 3. Typical nuclear industry devices 
such as resistance temperature detectors, and not so typical 
devices such as gamma densitometers, are used. The 
measuring devices numbered over 450 during the experiments 
examined in this study. Only the output of those which met 
the accuracy criteria of the testing directors are used in 


the comparisons made later. 


Table 3: Parameters Measured at the LOFT Facility 


- temperature 

- absolute pressure 

- differential pressure 
- material stress 

- liquid level 

- fluid density 

- flow rate 

- pump speed. 





Chapter 4 


MMS MODEL DEVELOPMENT 


Ame Development Process 


In order to create a model which suits the purposes of 
the systems analyst, it must first be decided which specific 
parameters he/she desires to examine. In this study, the 
parameters of interest include the temperature, pressure, 
and enthalpy at various points in the primary coolant 
system, and the temperature and pressure in the steam 
generator secondary side. These were chosen because the 
experiment reports included these parameters. Further, 
these parameters provide easy to interpret system 
performance factors which a person familiar with SReseniced 


water reactor operations will understand. 
4.2 Parameterization 


Once a module was selected for use in the LOFT facility 
model, the next step was to assemble input data which best 
describes the actual component. Collecting the necessary 
information required the use of many description documents 
and experiment reports. Each module used in the model 
required its own sources of information, and a set of 


calculations performed per the requirements of reference 
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[5]. The following sections describe the module selection 


process and the references used as sources for input data. 


4.3 Module Selection 


With the parameters of Section 4.1 in mind, module 
selection became a relatively (compared to the 
parameterization discussed later) simple process. The 
Criteria of selecting a module to represent an actual LOFT 
facility component now depended on what types of conditions 
a module was designed to simulate, and if those conditions 


were to be encountered in the transients used in this study. 


4.3.1 Major Components - Primary System 


The major components of the primary system include the 
reactor, steam generator, primary coolant pumps, and 
pressurizer. These were each selected based on the criteria 
presented in the module descriptions in Volume III of 
reference [5]. There are at least two modules available to 
represent each one of these components. In the transients 
used in this study, the LOFT facility's broken loop was not 
utilized, and so is not part of the models. Instead, a 


Single primary coolant loop was created. 
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Ameo. 1 Nuclear Reactor 


The current version of the MMS has four pressurized water 
reactor modules. The primary difference between the four is 
the number of core nodes used in solving the reactor 
kinetics equations. The simplest modules, RXl and RX1XR, 
use a Single node; RX3 uses three cylindrical stacked nodes; 
RX12 divides these three cylinders into four quadrants to 
provide twelve nodes. The RX1 module was designed for those 
transients which occur over periods of minutes and hours. 
Since decay heat plays a major role in the transients used 
in this study, it appeared that RX1l would not accurately 
predict the heat added to the primary coolant system after 
reactor shutdown. On the other hand, since most of the 
plant performance examined occurred immediately after a 
reactor shutdown, the more complicated set of twelve node 
kinetics equations was not considered necessary. Hence the 
module decided upon was RX3, a storage-resistive module. 
This module's kinetics were expected to be accurate enough 
for the short time the reactor was at power, and its decay 
heat calculations are the same as those of RX12. 

A feature of this module, natural circulation, also 
figured to be of use because in one of the transients 
examined the reactor coolant pumps were shut off. 

Input parameters for this module are the most complicated 
of any used in the LOFT facility model. Data from 


references [1], [3], [8], [9], and [12] were used to compile 
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the parameters and tables. Since the RX3 has two hot leg 
connections and four cold leg connections, identical to a 
Babcock and Wilcox two loop plant, one of the hot legs in 
the model was simply blocked with a boundary condition of 
zero mass flow. The cold legs, however, do not allow a zero 
mass flow rate because they input to the storage part of 
this module. Instead, the incoming cold primary flow was 
divided into four parts, one for each of the four RX3 cold 
legs. 

The name used for the reactor module in the LOFT facility 
model is RXX. Key internal variables include YRXX, ZATRXX, 
GBOR, and ZLHRXX, which are the rod heights, core power, 
boron concentration and upper plenum enthalpy, respectfully. 
Although the RX3 module allows for five rod banks (YRXX is a 
subscripted variable), at the LOFT facility the four actual 
control rods move as a Single bank. Hence all five values 
of YRXX were initially assigned the same value to simulate 


the single LOFT facility bank. 


cl 2) ee. Steam Generator 


There are five steam generator modules available in the 
MMS, one once-through version, and four U-tube versions. 
Since the LOFT facility steam generator is of the U-tube 
design, the once-though module was obviously not appropriate 


for use. 
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Of the four U-tube modules, one, UTSG, is a low order 
version. Low order versions are intended to be simple and 
relatively inaccurate to save on computation costs. Since 
the transients examined in this study involve large pressure 
changes on the secondary side of the steam generators, the 
low order module was originally not selected. Of the 
remaining three versions, two are intended to serve similar 
Functions: UTSGR and UTSGA. The latter 1S a newer, unproven 
module, while the former iS an unimproved version of the 
First U-tube steam generator module created. Improvements 
to UTSGR have been recommended by various users of the MMS, 
but have yet to be incorporated. The final version, UTSGE, 
includes feedwater preheaters that the LOFT facility does 
not have. 

For the reasons just explained, none of the U-tube steam 
generator modules were completely satisfactory for use in 
the LOFT facility model. Therefore two of the four modules 
avallable were selected to allow a wider range of MMS 
performance. One transient was examined using UTSGR to see 
if indeed the performance was aS poor as expected. UTSG was 
used for the second transient to determine if "low order" is 
a term that applies only during use in predicting extremely 
violent transients. 

These modules model the natural circulation of an actual 
steam generator, heat storage in the metal mass, and heat 


transfer by both subcooled and bulk boiling, but to varying 
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degrees. By its nature, the low order module uses a much 
Simpler set of equations. All of these characteristics were 
used to predict the performance of the LOFT facility under 
accident conditions. Both modules use the drift-flux method 
described in Chapter 2 to calculate the two-phase flow in 
their riser sections. 

UTSGR divides the U-tubes into four different heat 
transfer regions: hot and cold leg subcooled heating, and 
hot and cold leg bulk boiling heating. Different heat 
transfer equations are used, depending on the type of 
heating and temperature differences. The least amount of 
heat is transferred in the cold leg subcooled region 
because of the minimized temperature differences between 
primary and secondary fluids, while the most heat is 
transfered in the subcooled hot leg region because of the 
large temperature difference. UTSGR is a resistive-storage 
module on the primary side, and simply a storage node on the 
secondary side. 

UTSG uses a much simpler scheme for calculating the rate 
of heat transfer. Only two regions are used, one subcooled 
and one bulk boiling. Further, the sizes of these regions 
are fixed, while in the UTSGR the region sizes vary with the 
recirculation ratio. 

Both of these modules carried the named ITL for this 
Study. This name was selected because the only steam 


generator at the LOFT facility is in the intact loop. In 
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contrast to UTSGR, UTSG is a storage-resistive module on the 
primary side. Key parameters calculated by this module 
include the downcomer mass and mass flow rate. UTSGR also 


calculates a downcomer level. 


a a. 1. 3 Reactor Coolant Pumps 


The MMS has two pump modules, PUMP and PUMP4Q. PUMP can 
be powered by a variety of external sources including a 
Steam turbine and electric motor. Since the LOFT facility 
coolant pumps are driven by electric motors, this module 
would seem to be useful in this study. However, PUMP was 
designed to operate only in the positive differential 
pressure, positive flow region of the its operating 
characteristic curves. 

The PUMP4Q module does not have the option of an external 
power source. It does, however, simulate operation in all 
quadrants of the pump characteristic curves. Since one of 
the LOFT facility experiments used in this study included a 
shutdown and coastdown of the reactor coolant pumps, this 
module was selected. Another useful feature of this module 
is that it can Simulate two or more identical pumps in 
parallel without having to use a separate module for each. 
The name used for this module in the model is RCP, for 
reactor coolant pumps. 

Determining the input parameters of this module was 


difficult because not only were the volume and other 
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physical dimensions needed, but the actual pump's four 
quadrant operating curves were needed for input as tabulated 
data. Fortunately, this data was available in reference 


eZ). 


ao... 4 Pressurizer 


Of the four pressurizer modules available in the MMS, one 
could not be used because it is a low order version which 
does not continue operation when empty. The other three 
have similar characteristics, so the module settled upon was 
PZRB. This module 1S similar to one which was proven 
reliable in reference [6]. Implementation of this module 
waS Simply a matter of converting from the full-scale 
pressurizer physical parameters provided in the default to 
the much smaller dimensions of the LOFT facility 
pressurizer. The name used in the model for this module is 
PZ « 

The pressurizer modules are unique among those which 
carry fluids because they have three instead of two mass 
flow connections, one each for the surge line, spray line, 
and pressure relief line. Module operation ranges from 
nearly solid conditions to empty, a useful characteristic 
for studying reactor plant accident behavior. 

Key internal variables of PZRB are ZLSPZR, ZMLPZR, 
ZHTPZR, and ZWBPZR. These are, respectfully, the liquid 
level, liquid mass, electric heater power, and mass flow 


rate from the vapor to liquid regions. 





a 
as. 2 Balance of the Primary Coolant System 


The remainder of the primary coolant system includes two 
pipes, the surge connection, four valves, and two connective 


modules, depending on which steam generator module is used. 


ame 2. 1 Pipes 


To ensure proper loop flow time delays from the steam 
generator to the reactor and back to the steam generator, a 
Flow resistive pipe was placed between the reactor outlet 
and steam generator inlet (with a pass through the surge 
junction), and from the steam generator outlet to the 
reactor coolant pump suction. The latter pipe was not 
included with the UTSG steam generator because of the 
reversed location of the storage and resistive parts of this 
module compared with UTSGR. Instead, another hot leg pipe 
located between the surge connection and the steam generator 
primary coolant inlet was required. This module is named 
SSG, for surge to steam generator connection. When UTSG is 
used, the volume of the hot leg pipes iS increased as 
required to account for the entire volume of all piping. 
Proper accounting for the correct volumes is expected to 
ensure proper pressure calculations. 


The hot leg pipe is called RXO for reactor outlet. 
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fees 2. 2 Sumge dune ton 


The surge junction acts as a storage node to the flow 
From pipe RXO into the primary side of the steam generator, 
and as a resistive node to the flow into and out of the 
Surge connection of the pressurizer. Although very little 
energy is lost or gained in this node, a separate surge 
module was developed for the MMS and is used in the LOFT 
Facility model to account for the large reactor coolant 
Flow, much smaller surge flow, and to provide a resistance 
to keep the presSurizer from completely emptying into the 
rest of the primary loop during an up-power transient. 
Normally, at least two separate modules would be required to 


allow bi-directional flow. 


al, BE Valves 


The four valves used in the primary coolant model are the 
pressurizer spray valve, PSY, the pressurizer safety relief 
valve, REL, the high pressure injection stop valve, XC, and 
a valve, BRK, which represents a connection between the cold 
leg and the blowdown suppression tank. 

There iS a pipe, with a gate valve, on the actual LOFT 
Facility cold leg piping downstream of the reactor coolant 
pumps. It is connected at its other end to the blowdown 
Suppression tank. Flow through this connection is required 


in one of the transients examined in this study. In the 
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LOFT model only the valve is required. The downstream side 
of the valve has simply a boundary condition, a constant 
pressure which was varied during model execution from 14.7 
to 500.0 psia. The need for BRK will become apparent in 
Chapter 5. 

The spray control valve PSY connects the primary coolant 
pump discharge to the vapor space of the pressurizer. It is 
an automatically controlled, quick opening valve both at the 
LOFT facility and in the MMS model. 

REL connects the presSurizer vapor space to the blowdown 
Suppression tank at the LOFT, but simply discharges to an 
infinitesimally large tank in the LOFT facility model. The 
pressure in this "tank" iS a constant 14.7 psia by use of a 
boundary condition. Verification of the use of this number 
was not made because in the transients of this study, no 
Significant pressure increases occurred. This fact allowed 
the model to be completed without modules representing the 
blowdown suppression tank and associated piping. REL, like 
PSY, 1S modeled as an automatically controlled and quick 
opening valve. 

Finally, XC is also an automatic quick opening valve. It 
is the only part of the Emergency Core Coolant System used 
in the model. Since the high pressure injection flow is 
provided by a positive displacement pump, the actual 
upstream pressure of this valve always follows the 


downstream pressure. This effect was included by making the 
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pressure upstream a constant 20 psid greater than the cold 


leg pressure. 


Ames 2. 4 Connections 


Three other modules, a flow divider called RXI, a 
Junction called JUN, and a simple connection called PRX, 
complete the primary loop. RXI, located at the inlet of the 
reactor, divides the inlet flow into six streams: four are 
the reactor cold legs described in Section 4.3.1.1, one is 
the pressurizer spray line, and the sixth connects to the 
break connection valve, BRK, described earlier. 

The junction module, JUNC, is required for the high 
pressure injection system connection. At the LOFT facility 
this connection, along with that of the pipe to the blowdown 
Suppression tank, is at the coolant pump discharge. In the 
MMS model, two connections are needed because the injection 
Flow is into the primary coolant system, while the flow to 
the blowdown suppression tank is out of the system. The MMS 
does not allow flow into and out of a system boundary at the 
Same point. 

The connection module, PRX, is needed because the pump, 
RCP, is a storage-resistive module, and would otherwise 
discharge directly into the reactor, a resistive-storage 
module, leaving two adjacent resistive nodes. (The dividing 
node is neither storage nor resistive). The MMS 


incompressible fluid connection module, appropriately 
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invoked by the command CONNI, is a storage module with a 
volume of Zero. 

The 14 modules of Section 4.3 which make up the primary 
coolant system, their joint configuration, MMS module names 
and LOFT facility model names are shown in Figure 14. The 
steam generator module of this figure is UTSGR. Figure 15 


Shows the substitutions made to use UTSG. 


Ae 3. 3 Secondary Coolant System 


The secondary coolant system is modeled by up to five 
modules, in addition to the secondary side of the steam 
generator. These modules include three valves, one flow 


divider, and a connector. 


Amesis 5 5 Steam Control Model 


The LOFT facility steam flow rate was varied by manual 
and automatic control in the transients examined in this 
study. The MMS, however, has no provisions for both types 
of control on the same valve. Therefore, two main steam 
control valves were used in the model, although the LOFT 
facility has only one. Under normal and anticipated 
transient conditions, both were not expected to be open at 
the same time, so that an abnormally high steam flow rate 
would not exist. Identical valve modules that allow for 


compressible flow, VALVEC, are used. The manually operated 





VALVEI- 
“BRK” 


Figure 14: 
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Figure 15: LOFT Facility Model (Substitutions made for 
UTSG) 


valve is called MSS for main steam stop, and the 
automatically operated valve is MSR for main steam relief. 
At the LOFT facility, the relief function is accomplished by 
automatic operation of the normal flow control valve. 

Flow from the steam generator to the steam valves 1s 
through the divider MSL to allow flow directly to whichever 


of the valves 1S open. 


3, oe ry Feed System Model 


The feed system consists of only two modules, the feed 
control valve FRV, and an incompressible connection RFW. 


Feed flow rate is normally controlled by the position of the 
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Figure 16: Secondary Coolant System Model 


feed control valve. However, to accurately recreate the 
feed flow rate at the LOFT facility, RFW was added. To 
determine the effects of valve position on flow rate would 
have required an extensive testing procedure. Instead, use 
of a connective module allows flow rates to be input as a 
boundary condition. 

The secondary coolant system model and module names are 
shown in Figure 16. Appendix B is a listing of all the non- 
control module names and their equivalent LOFT facility 


components. 
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ees. 4 Control Components 


The only controllers needed in the LOFT facility model 
are on-off switches and their associated actuators. All 
automatically controlled valves require both a switch and an 
actuator module. The pressSurizer heaters and the low 
pressure sensing reactor shutdown switch require only an on- 


off module. 


aes. 41 Pressurizer Heater Control 


The LOFT facility pressurizer has two sets of immersion 
type electrical heaters for automatic pressure control, the 
12 kw cycling heaters and the 36 kw backup heaters. Each 
cycles on and off controlled by the pressurizer pressure. 
In the model this pressure is the variable PPZR. The MMS 
pressurizer modules allow for direct control of the heaters 


based on any system pressure selected by the modeler. 


woes. 4.2 Rod Position Control 


The low pressure automatic reactor shutdown or scram 
Switch operates Similar to the heater controllers, but with 
rod position as their output. This is the only automatic 
shutdown used in the model because it is the only one which 
occurred at the LOFT facility during the selected 


transients. 
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fea. 3 Valve Control 


Table 4 lists the valves in the LOFT facility model which 
operate automatically, the parameter which controls them, 
and their associated actuator. These valves are part of 
both the primary and secondary systems, and perform a 
variety of functions. The other model valves are controlled 
by boundary value tables. 

The total number of modules used in the LOFT model is 30: 
14 in the primary coolant system, five in the secondary 
coolant system, and 11 control modules. Appendix C is the 


fully assembled LOFT facility model in MMS format. 


4.4 Assembly Process and Initialization 


Once the modules were selected and the initial parameters 
set, each was "operated" individually using as input the 
initial conditions of one of the transients described in 
Chapter 5. When a module operated satisfactorily in the 
correct steady state by itself, it was set aside and the 
next module was tested. 

Assembly of the complete model began with the reactor, 
RXX, operating alone at 100% power, flow rate, enthalpies, 
and pressures. Next the hot leg pipe module, RXO, was 
attached. The now two-module model was operated until it 
worked satisfactorily in the steady state, and then another 


module was attached. This process of attaching a module and 
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Table 4: Model Valve Controllers 


Valve On-off Module CeneHOLling Associated 
Name Parameter Actuator 
Spray valve SVC PPZR PSY 
PSY 
High Pressure XCC PPZR XC 
Injection XC 
Pressurizer RVC PPZR REL 
Relief REL 
Main Steam MSC Steam Dome MSR 


Relief REL Pressure PSTO 
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then testing to ensure the proper steady state parameters 
were calculated was repeated until the primary loop, without 
the connection from pump discharge to the pressurizer, was 
completed. The last module attached to complete the loop 
was the steam generator. Testing with this module in the 
system required that some of the constants, in particular 
the temperatures of the U-tube metal, be altered slightly 
From their calculated values. Since the calculated values 
of these temperatures were educated guesses at best, 
changing them slightly was not treated as a significant 
problem. When completed, the loop's heat balance, flow 
rates, and pressures were correct. 

Next the spray valve connection was made to the reactor 
inlet piping. With this connection in place, the model 
would not operate at all. After a lengthy investigation 
into the problem, it was discovered that if a small flow 
resistance due to shear stresses was input into a relatively 
small diameter pipe, the model would stop execution. [In 
this case, the resistance of the spray line was much less 
than that of the much larger sized reactor coolant inlet 
piping. The stop occurred because all flow would attempt to 
go through the much smaller spray line, causing flow 
reversals in the large dimension reactor inlet pipes. The 
MMS terminates execution if it senses such flow reversals. 
This problem was solved by decreasing the flow conductance 


(and thereby increasing the resistance) term of all the 
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small diameter pipes and valves in the model. RXO and SGP, 
the pipe modules representing the much larger hot and cold 
legs, respectfully, did not require large changes in their 
flow conductances. Only minor adjustments brought their 
differential pressures to the required values. 

With the primary coolant carrying modules now all 
attached, the next step was to add the control modules. 
This step, also, was accomplished by adding one module at a 
time, followed by testing to ensure the entire model still 
operated satisfactorily. Addition of the pressurizer heater 
and spray controllers caused a somewhat awkward problem in 
that both the transients examined began with pressurizer 
pressure at about 2158 psia. The LOFT facility controllers, 
and those in the MMS model, are set to control pressure at 
2250 psia. As Figure 17 shows, allowing the model to reach 
steady state causes some short fluctuations followed by a 
steady rise to the design pressure. Figure 18 shows the 
action of the heaters during the period that the model 
searches for its steady state. Note that all heaters are on 
until 350 seconds, when the backup heaters are de-energized. 
At steady state, the model runs continuously with pressure 
just below the turn-off set point of the cycling heaters. 

Finally, the secondary modules were attached, using the 
same add and test process. The only significant adjustment 
required to the previously calculated secondary parameters 


was to the flow conductance of the main steam control valve, 
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MSS. This adjustment was needed to ensure that the steam 
flow rate out of the generator was maintained at 808 psia. 
The constant parameters which required changes are shown 
as the first series of lines in the ACSL command files of 
Appendix D. This listing illustrates that to change an 
input parameter's value that has already been set by a 
FORTRAN statement, all that is required is a SET statement 
in the command file. Modification of the ACSL program, 


retranslation, and recompilation are not required. 





Chapter 5 


TRANSIENT DESCRIPTIONS 


Slt Selection Criteria 


In order to validate operation of the MMS model of the 
LOFT facility under a variety of conditions, two transients 
were selected for analysis. The basic requirements to be 
met for inclusion of a transient in this study are listed in 


Table 5. 


Table 5: Transient Selection Criteria 

1) Adequate data had to be available for the key 
thermodynamic properties of both the primary and 
secondary coolant systems to at least 200 seconds 
after transient initiation; 

2) The transients should include primary system 
pressure or power changes that initiate high 
pressure injection and/or a reactor scram; 

3) Comparison data to another thermal-hydraulic 
reactor analysis code should be available. 


Several of the LOFT facility experiments meet two or more 
of these criteria, so the selection was further refined to 
those transients which were initiated from easily achievable 
Steady states. These steady states are those with the 
reactor at a constant power in the power range. Transient 
initiation from below the power range would have required 


either a steady state achieved from a previous scram 
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transient or elimination of the RX3 module, and use of an 
appropriate substitute as a source of decay heat. A steady 
State achieved from a previous scram transient was not 
acceptable because the ability of the MMS to model a LOFT 
facility scram had not been validated. As described in 
Chapter 1, this validation is the objective of this study. 
Two experiments were chosen for analysis: 4L3-5, a small 
break loss of coolant experiment, and L6-3, a rapid rise in 


secondary coolant capability. 


SZ Transient Procedures and Significant Events 


pee. L Small Break Loss of Coolant Experiment, L3-5 


L3-5 is one of a series of six small break experiments 
performed at the LOFT facility soon after the Three Mile 
Island accident of 29 March 1979. The objectives of this 
series include: 


To determine the important plant thermal, 
hydraulic, operational, and neutronic phenomena 
during a variety of small break LOCEs (Loss of 
Coolant Events). . . 

To evaluate the effectiveness Gmrnews “Ss in 
mitigating a slow 

depressurization transient. ; 

To determine the effect of primary coolant pump 
operation_on plant 

response. 


7 Leanne Thuy Lien Dao and Janice M. Carpenter, Experiment 
Data Report For LOFT Nuclear Small Break Experiment 
L3-5/L3-5A (Idaho Falls, Idaho: EG&G Idaho, Inc. 1980), 


jane Sie 
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The objectives of the L3-5 phase were: 

1. To conduct a small break depressurization 
in the LOFT facility with a 16.19-mm (0.6374-in.) 
diameter break orifice in the intact loop cold leg 
between the primary coolant pump and the reactor 
vessel, with primary coolant pump trip at the 
rupture, with the HPIS injecting into the reactor 
vessel downcomer, and with the accumulator 
isolated from the intact loop 

Ze. To measure the primary system coolant 
inventory and system coolant mass distribution as 


a function of time during the depressurization 
uSing available instrumentation. 


5.2 [at Initial Conditions 


The initial operating conditions at the facility for 
experiment L3-5 are presented in Table 6. These data are 
the actual measured values of the parameters listed, which 
are not necessarily those values intended by the operators. 
However, all the initial conditions presented are within 2 
of their specified values, except for pressurizer liquid 
level and steam generator water level. These were 0.12 m 
above and 0.06 m below their specified values, respectfully. 
Listed with the actual values in Table 6 are those 
calculated by the MMS model after it achieved steady state. 
There are Significant differences between the actual and MMS 
values of steam generator level and boron concentration. 
The boron concentration difference was ignored because the 
reactor was shutdown by a scram 4.8 seconds before transient 


initiation, and only decaying fission products affected the 


8 Dao and Carpenter, pp. 3-4. 





Table 6: Initial Conditions for L3-5 


Parameter Initial Values 
LOFT MMS 
Primary mass flow rate 3.77 x 106 3.78 x 10° 
(lbm/hr ) 
Hot leg temperature 577.4 577.8 
le) 
Cold leg temperature 545 5427 
ee) 
Pressurizer pressure ZUG 2 7a AIBN 
(psia) 
Pressurizer level 4.16 4.18 
Ore), 
Boron concentration 650 1350 
(ppm) 
Control rod height 83 81 
(% withdrawn) 
Reactor power level 49 50.8 
( MW ) 
Secondary mass flow rate 209088 224730 
(lbm/nhr ) 
Steam generator level 10.3 16.83 
Ob), 
Steam generator pressure 809 808 


(psia) 
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heat generation within the reactor. However, the steam 
generator level problem may have had some effect on the 
MMS's ability to predict heat removal performance, since the 
higher downcomer level indicates a larger liquid mass on the 
secondary side. 

Attempts to change the liquid level by changing the 
physical parameters used in the input calculations only 
succeeded in radically altering other steam generator 
parameters, including the rate of heat transfer out of the 
primary system. Since the heat balance on the steam 
generator with this water level was very nearly correct, it 
was decided to continue the transient calculations from this 
pormet. 

Since the highest rates of change of the parameters of 
interest occurred within the first 200 seconds of this 
experiment, this is the only period examined in this study. 
However, reference [4] contains experimental data to 2400 


seconds if needed for future investigations. 


ee2sl.2 Significant Events 


Preparations for experiment L3-5 began by taking the 
reactor critical about 45 hours prior to transient 
initiation, and raising the power level to 49.3+.7 MW at 
35.9 hours prior to initiation. This power level was 
maintained to allow a near equilibrium buildup of decay heat 


generating fission products. Such an equilibrium would then 
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Simulate a power station operating many days at near 100% 
rated power. 

At seven minutes prior to initiation, the instrumentation 
and recording systems were turned on, and at T=-4.8 seconds, 
the reactor was manually scrammed. With the indication of 
four “rod bottomed" lights at T=-2.8 seconds, the experiment 
was begun. (All times are referenced to T=0, the point at 
which the leak was begun.) 

At T=0 seconds, the small leak was simulated by opening a 
valve in the drain line between the reactor coolant pump 
discharge in the cold leg and the blowdown suppression tank. 
This is the valve modeled by BRK in the MMS model. The 
initial flow rate into the blowdown suppression tank was 
43000 l1bm/hr. This rate decreased rapidly to 25000 
lbm/sec. at T=150 sec. At this point the flow rate out of 
the primary system slowly dropped until the experiment was 
concluded at 2309 sec. after initiation. The leak was sized 
to simulate a four-inch-diameter break at an actual 3000 
MW(t) PWR plant. Sizing was accomplished by using an 
Orifice in the drain line. 

At T=0.8 seconds, the reactor coolant pumps were tripped. 
The pumps coasted down until their motor breakers tripped 
open at T=17.7 seconds. This marked the end of the coast 
down and the end of forced circulation of the primary 
coolant. Interestingly, natural circulation was detected as 
soon as T=17 seconds, indicating that flow never fully 


stopped. 
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At T=4.0 sec., HPIS flow was initiated into the reactor 
downcomer, and continued until past the time of interest in 
this study. In the MMS model of the LOFT facility, this 
Flow is into the reactor coolant pump discharge. Since the 
ultimate effect of the HPIS flow is to mitigate pressure 
losses and to maintain coolant inventory, the location of 
injection should not have made a difference in the MMS 
predictions. Automatic initiation of injection began when 
the primary system pressure reached 1915 psia. None of the 
references used in this study indicated precisely where this 
pressure was measured, so it was assumed to be in the intact 
loop hot leg. In the first four seconds of this experiment, 
pressure dropped 243 psia to reach the HPIS injection set 
point. Such a rapid change was caused by both the coolant 
Shrinkage due to the scram and the loss of coolant through 
the simulated leak. A large amount of heat was still being 
drawn off by the steam generator because the main steam 
isolation valve requires about 10 seconds to shut. This 
valve began to close shortly after the scram, but was not 
fully closed until T=5.2 seconds. 

At T=22.2 seconds, the pressurizer was emptied. Pressure 
in the pressurizer at this point was 1450 psia. About 6.2 
seconds later, the reactor's upper plenum reached saturation 
conditions of 572°F and 1250 psia. Here the pressure drop 
Slowed considerably. At T=30.0 seconds, hot leg voiding 


began, and at T=80.0 seconds, cold leg voiding began. 
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Finally, at T=92.9 seconds, flow through the leak reached 
Saturation. At T=200 seconds the primary system pressure 
was down to 986 psia. 

The only other significant event to occur in the first 
200 seconds of L3-5 was the automatic initiation of 
auxiliary feed flow at T=63 seconds. This flow continued 
until well past the 200 second mark. In the MMS model, this 
flow is provided by a table of boundary values. 

Figures 19, 20, and 21 show the changes of some 
Significant parameters of this experiment, along with the 


changes predicted by the RELAP code, if available. 


Gaeeee, 'xcessive Steam Load Experiment, L673 


The objectives of the L6 series of experiments include: 


- determine the important thermal, hydraulic, 
Operational, and neutronic phenomena during an 
anticipated transient at the LOFT facility and to 
identify any unexpected behavior... 

- provide data to evaluate reactor transient 
analysis techniques used to analyze anticipated 
transients .. . 

- provide data to assist in analyzing the 
relationship between behavior in LOFT and in, a 
commercial PWR during anticipated transients. 


The specific objectives of the excessive steam load 
experiment, L6-3, were: 
a. Investigate plant response to a transient in 
which the heat removal capability of the secondary 
system is significantly increased 


b. Provide continued evaluation of automatic 
recovery methods 


9 Nalezny, p. 9-l. 
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Figure 19: Experiment L3-5: Hot Leg Flow Rate 
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Figure 20: Experiment L3-5: Primary Coolant System Pressure 
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Figure 21: Experiment L3-5: Pressurizer Level 
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c. Provide data to evaluate code capabilities for 
secondary system 
initiated events.19 


pee, 2. Initial Conditions 


The initial conditions for this experiment were much the 
same as those of L3-5, with the major exception being that 
the reactor was at 75% of rated power instead of at maximum 
power. The actual initial conditions are listed with those 
of the MMS model for this experiment in Table 7. 

Of significance is the lower steam generator pressure, 
775 psia at 75% power vs. 809 psia at 100% power. This 
lower pressure was difficult to recreate in the MMS model 
because reducing the steam flow by closing the main steam 
control valve caused an increase in the model's steam 
generator secondary side pressure. This is, of course, what 
would happen in an actual reactor plant without some sort of 
automatic primary temperature control. However, the average 
Primary temperature in the U-tubes was 548.5°F at 75% 
power, compared to 560°F at 100% power. This indicated 
that the facility's reactor has a load following primary 
temperature control system, although such a system was not 
described in any available reference. Modeling this system 
was not required to meet the objective of this study, but 
the changes from the L3-5 initial conditions to the L6-3 


initial conditions were required. Hence, an entirely new 


10 Walezny, Dp. 9-2. 





Table 7: Initial Conditions for L6-3 


Parameter 


Primary mass flow rate 
(lbm/hr) 


Pressurizer pressure 
(psia) 


Cold leg temperature 
hot ) 


Hot leg temperature 
ee) 


Reactor power level 
( MW ) 


Control rod position 
(% withdrawn) 


Pressurizer level 
Et ) 


Secondary mass flow rate 
(lbm/hr ) 


steam generator pressure 
(psia) 


Initial Condition 


LOFT 
3.80 x 10° 
PUES 6 
535.3 
561.7 
36.9 

81 

3.94 
163944 


tiie, 


MMS 





3.85 x 10° 
2194.7 
538.2 
564.8 
38.0 

81 

3.90 
163940 


UG S56 


81 





82 


set of steady state parameters was developed. The model's 
primary temperature was lowered by raising the boron 
concentration to 1397 ppm, or 52 ppm more than the level for 


L3-5 e 


ez a2 Significant Events 


Criticality for experiment L6-3 was achieved about 16 
hours prior to experiment initiation. At about 4.5 hours 
prior to initiation, power was raised to 49.5 MW, and then 
lowered to 36.9 MW (75% rated power) just before the 
experiment was begun. 

L6-3 was initiated at T=0 seconds by ramping open the 
steam flow control valve from the 75% power position. As 
steam flow increased the cold leg temperature dropped, 
cauSing an increase in the reactor power level due to the 
negative temperature coefficient of reactivity. Power 
reached a maximum of 42.2 MW at 15.6 seconds, when the 
reactor scrammed. Scram was initiated automatically upon 
receipt of a low primary system pressure signal at 2080 psia 
in the intact loop hot leg. 

The pressure decrease continued because of delays in 
reducing the steam flow. The steam flow control valve was 
not completely shut until T=36.2 seconds. At the LOFT 
Facility, this valve shuts automatically upon receipt of a 
reactor shutdown signal. Delays in reaching the fully shut 


position include time for a "close" signal to be sent to the 
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steam flow control valve actuator, time to reverse the 
direction of travel of the control valve, and the travel 
time to reach the fully shut position. These delays totaled 
about 20 seconds. The feed pump was tripped immediately 
after the scram. 

The HPIS pumps started automatically at T=26.4 and T=26.6 
seconds when primary system pressure reached 1915 psia. The 
pressure drop was immediately mitigated, and within ten 
seconds, pressure began to rise again. The HPIS pumps were 
shut off by the operators at T=48.6 and T=50.0 seconds, with 
pressure at 2100 psSia. By the end of the period of 
interest, 200 seconds after experiment initiation, pressure 
had almost returned to the automatic control band. 

Decay heat input to the primary coolant was near its 
maximum possible rate because of the previous operating 
history near 100 power. The reactor decay heat generation 
exceeded the steam generator heat removal at T=33 seconds. 
This effect, too, helped to quickly restore pressure to its 
normal level. 

Figures 22, 23, and 24 show the changes in steam demand, 
reactor power, and primary system pressure, respectfully, 
that occurred during this experiment. Also in these figures 


are the applicable RETRAN predictions. 
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Figure 22: Experiment L6-3: Steam Demand 
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Figure 23: Experiment L6-3: Reactor Power 


5.3 Applicability to the Validation Process 


It should be noted that both of these transients involved 
large changes in primary pressure after a reactor scram. 
However, experiment L6-3 was not as "violent" as L3-5 
because no loss of coolant occurred, the primary coolant 
pumps were not shut off, and saturation conditions did not 
occur in the hot leg piping. These two experiments, then, 
compliment each other in the MMS validation process. It was 
considered possible that the MMS could handle predicting the 
LOFT facility performance in L6-3, yet not be appropriate 


for such severe transients as those of L3-5. 
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Figure 24: 


Experiment L6-3: Primary Coolant System Pressure 





Chapter 6 


PERFORMANCE OF THE MODULAR MODELING SYSTEM 


OL Experiment Predictions 


Operation of the MMS models produced output consisting of 
many modeler selected thermodynamic state variables and 
module internal variables. Many executions were required to 
trim the model input so that reasonable output could be 
obtained. The first sections of this chapter describe the 
best model performances in predicting the parameters of 
experiments L3-5 and L6-3. The last section is included to 
demonstrate some of the problems that can be encountered 


when working with a code such as the MMS. 


oo. 1 Experiment L3-5 


The MMS predictions of the LOFT facility's thermal- 
hydraulic performance in experiment L3-5 begins with the 
reactor plant at steady state. The steady state operating 
parameters are listed in Table 6. To execute this 
transient, tables of various operator actions were made part 
of the ACSL program. The variables changed by these tables 
are the steam flow control valve position, the feed flow 
rate, and the position of BRK, the leak simulation valve. 
The L3-5 time references in this chapter begin at 10 seconds 


before the opening of the "leak" valve, BRK. Hence the 
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scram occurs at T=5.2 seconds, and BRK opens at T=10.0 
seconds. 

The parameters considered key in evaluating the transient 
performance of the MMS model are the pressurizer level, the 
primary system pressure, the primary coolant flow rate, and 
the steam generator secondary side pressure. These were 
selected because they can be used both as direct performance 
indicators and have the synergism to be used for 
interpreting other parameters. These other parameters 
include temperatures and mass inventory. Figures 25 through 
28 show the actual reactor plant's trends of the four 
selected parameters compared with the MMS predicted values. 
The actual values are indicated by the symbol "{4j," the 
predicted values by the "«," 

As these figures readily indicate, the MMS did not 
Operate past 30 seconds into the transient. Discounting the 
first ten seconds used to adjust the zero marker and allow 
for the scram, less than 20 seconds of the actual transient 
are shown. Of approximately 200 attempted reinitiations and 
executions of experiment L3-5, the longest "real" time that 
was reached was achieved on the runs shown in Figures 25 
through 28. It should be noted that a typical MMS run uses 
about 1 CPU second for each second of real time up to about 
20 seconds, and 1 CPU second for 10 real seconds thereafter. 

As shown in Figure 25, the MMS appears to accurately 


predict the changes in the primary coolant cold leg flow 
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rate. By the built-in variable naming convention of the 
MMS, WCPO indicates this is a plot of the mass flow rate, W, 
at the location named by the modeler as CPO. In the LOFT 
Facility model, CPO is the discharge of the primary coolant 
pumps. After the coolant pumps are shut off, pump coastdown 
takes a number of seconds, during which flow steadily 
decreases. The largest difference in flow rates occurs at 
T=28.0 seconds. Note that because the pump speed was 
entered as a series of steps using the ACTION command in the 
ACSL command file, the flow coastdown took a scalloped 
appearance. The ACTION command allows variables to be 
changed more than once in a Single model run. However, the 
changes are instantaneous, as opposed to the smooth 
appearance of changes made using a TABLE command. The TABLE 
command, on other hand, requires recompilation of the entire 
model if even a single value is changed. Use of the ACTION 
command allowed frequent changes, at minimal cost, in the 
pump coastdown rate in order to extend the the model's 
Operating time. The coastdown in Figure 25 is at the actual 
rate of the facility's coolant pumps. Changing the 
coastdown rate did not change the time at which model 
execution terminated. 

Figures 26 and 27 show together the changing of the 
pressure in the primary system, with the accompanying 


lowering of the pressurizer water inventory. 





4.80 


Abe AR MIP id 


h Au 


meoure 26: 


L3-5 Predicted Performance: 
Level 


Pressurizer Water 


ou 





a2 


2.20 ZOU 


Ik gsi 


PCPI «103 (PSI) 
1,60 


aU 


1. 00 


0.00 0) 76 f. at lea) a0) 8 


a0 30.0 
ZIV IME (oie) 


megure 27: L3-5 Predicted” Performance: Primary System 
Pressure 





a3 


ZLLPZR and PPZR are the MMS names for pressurizer level and 
pressure, respectfully. The predicted level and pressure 
drop slowly after the scram, and at an increasing rate when 
the leak is initiated. The drop in the level follows the 
actual rate initially, but slowly decreases due to the high 
pressure injection flow until the levels are equal at T=25 
seconds. High pressure injection was also initiated at the 
LOFT facility, of course, but its effects seemed to be 
masked by instrument inaccuracies. Execution termination 
occurs shortly thereafter, with an interesting sharp upturn 
in the predicted level. This upturn is the leading clue to 
the cause of the failure of the model to continue past 30 
seconds. This failure will be discussed further in Chapter 
7. The pressurizer pressure diverges from the actual value 
at the time of leak initiation. The predicted pressure 
change does not reach the high rate of change measured at 
the actual reactor. At T=12.5 seconds, shortly after 
initiation of high pressure injection, the rate of change of 
actual pressure slowed. Here the predicted and actual rates 
became nearly equal. The error between predicted and actual 
values at 28.5 seconds is 11.7%. 

Finally, Figure 28 shows the pressure in the steam 
generator begins to increase after the flow control valve is 
shut, as expected. The model predicts a much faster initial 
rise in the pressure, but then appears to slowly approach 


steady state. Of course, what happens after 30 seconds is 
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not known, but it appears that the predicted steam generator 
pressure would never reach the relief valve setpoint, while 


plant data shows that the actual relief valve did open. 


Ore)... 2 Experiment L6-3 


Figures 29 through 31 show the trends of the key 
parameters of experiment L6-3. For this less severe 
transient, the primary flow rate is not shown because the 
pumps were not shut off, nor is the steam generator 
secondary side pressure. Changes in this pressure were 
instead used as input data, along with the secondary flow 
rate, to initiate the transient on the model. Added to the 
figures is reactor power, shown in Figure 29. 

Since this transient began with the opening of the steam 
flow control valve at a time of 0 seconds, the figures used 
for the evaluation also begin with a time of T=0 seconds. 

Of immediate note is that the MMS operated to the 
intended end of the transient, 200 seconds. The power trend 
is plotted for only the first 50 seconds because once the 
reactor was shut down, the MMS predicted power remained 
consistently about 5% above the the actual plant data. The 
50 second plot expands the first part of the transient for 
better clarity. 

The key differences between the actual and predicted 
performance of the reactor power are the rate of change of 


power while steam flow was increasing and the maximum power 
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Figure 29: L6-3 Predicted Preformance: Reactor Power Level 
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reached before the reactor is automatically shut down. The 
MMS rate is higher and thus when the reactor scrams on low 

primary system pressure the model power peaks at 100% while 
the actual power was 86%. If the model had also included a 
high power scram setpoint, it may have been reached before 

the reactor scrammed on low system pressure. 

In Figure 30, it can be seen that the pressurizer 
pressures initially track very closely together through the 
time of the scram and up to about T=20 seconds. Both the 
LOFT facility and MMS low pressure scram set points are set 
to initiate reactor shutdown if pressure drops past 2130 
psia in the hot leg. In experiment L6-3, as in experiment 
L3-5, the predicted pressure drops more slowly than the 
actual value. 

The pressure drop continues until the initiation of high 
pressure injection. Because of the slower pressure drop 
predicted by the MMS, HPI does not come on until T=36 
seconds, 10 seconds after the actual initiation. The slow 
reversal in the rate of change from negative to positive is 
predicted to occur more rapidly than did the actual 
reversal. Hence the minimum pressure reached by the model 
is 1960 psia, while in reality, pressure reached a minimum 
of 1940 psia. Of note is the almost immediate mitigation of 
the rapid pressure change by the model. 

With the HPI system running, the predicted pressure rises 


at nearly the rate of the actual increase. The actual HPI 
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pumps were turned off by the operators, causing the sharp 
change in the rate of pressure increase when pressure had 
returned to 2100 psia. This pressure is 70 psia below the 
automatic pump shutoff set point. Using only an on-off type 
of controller for the HPIS, both automatic and manual pump 
operation was not allowed. Hence the model's pumps are not 
turned off until primary system pressure returns to 2170 
psia. At this point the predicted pressure rise slowed, and 
primary system pressure control was eventually returned to 
normal operation. At T=200 seconds, the predicted pressure 
was back in the normal operating range, while the actual 
pressure had reached only 2105 psia. 

The predicted water level of the pressurizer tracked very 
closely with that of the actual experiment. Again, because 
of the slower drop in the predicted pressure, the time at 
which the level is predicted to reach itS minimum is a few 
seconds after the actual time. The final predicted level is 
above that of the actual level with an error of 7.8%. 
Figure 31 shows the actual and predicted values of the water 


level. 


6.2 Unsatisfactory Model Execution 


As described earlier in this chapter, many executions of 
the MMS models were required to arrive at the results of 


mtiamres 25 through 31. Figures 32 through 35 show the 
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results of some unsatisfactory runs. Although these 
executions are termed unsatisfactory, they were not without 
use. Most provided some indication of a cause of execution 
failure or a required adjustment in one of the model's input 
parameters. 

Figure 32 shows the change in pressurizer level on an 
early attempt at predicting experiment L3-5. The pressure 
shown is that at the primary coolant pump suction. This 
pressure is initially, as expected, about 20 psia less than 
the presSurizer pressure. 

After the scram at T=5.2 seconds, the pressure drops a small 
amount, but then begins a slow climb past its initial value. 
At T=35 seconds, the pressure again began to drop, but at a 
Slowly decreasing rate. Investigation showed that the 
ACTION command used to initiate the scram was inserting only 
one of the five simulated control rods banks, while the rest 
remained at their critical position of 53.5 inches. The 
model's reactor became momentarily subcritical, causing more 
primary system heat to be removed by the steam generator 
than was being input by the the reactor. This heat removal 
difference caused the small drop in primary system pressure 
between T=7 seconds and T=13 seconds. Because the primary 
coolant's average temperature dropped, the reactor returned 
to a critical condition due to the negative temperature 
coefficient of reactivity. At T=13 seconds, the heat 


removal rates were nearly equal, and pressure began to 
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increase because the backup heaters came on. Finally, when 
the backup heaters deenergized, the pressure began to fall 
quickly due to the coolant leak, causing the heaters to come 
on once again. This run demonstrated the unSuitability of 
uSing the ACTION command in changing the value of a 
Subscripted variable. The rod heights of module RXX are 
contained in the five values of YRXX. The ACTION command 
changed only the value of YRXX(1). In order to change the 
remaining four values of this variable, the actual ACSL 
generated FORTRAN code had to be modified. 

Figure 33 shows the steam generator secondary pressure of 
the facility, uSing the initial conditions of experiment 
L3-5. In this simulated run, however, the leak was not 
initiated, nor were the reactor coolant pumps secured. The 
pressure in the steam generator increases at nearly the same 
rate as the actual pressure, until the steam relief valve 
opened. The MMS valve was modeled as quick opening, since 
no other description was found in the references. This plot 
indicates that the actual valve has some level of 
accumulation, and does not fully relieve the pressure until 
it has reached a more fully open position. Hence the 
predicted pressure drops quickly upon opening the relief 
valve, and then builds up again when the valve shuts. Two 
valve operation cycles occur before the steam generator heat 
removal has been reduced to a point where the secondary 


pressure no longer reaches the valve's opening setpoint. 
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The final trends are similar to the actual changes in the 
LOFT facility steam generator pressure. 

A number of L3-5 runs were made with varying primary 
system volumes and varying leak flow rates. The rate at 
which the coolant flowed through BRK had a direct impact on 
how fast the primary pressure dropped, as expected. It was 
also expected that changing the system volume would change 
the rate of depressurization. However, this was not the 
case. Even when the volume of the hot and cold legs was 
reduced to less than a tenth of the actual value, the rate 
of change of the pressure was not affected. The reason for 
this effect is unknown. 

The effects of varying the flow rate of the high pressure 
injection pumps are shown in Figure 34. Here the initiation 
set point pressure was too low, moving the predicted 
pressure increase curve to the right of the actual curve. 
Once the HPI pumps did come on, however, they very quickly 
returned pressure to above the inital value, up to the 
automatic control band. This plot showed that the model's 
HPI flow rate was too high, and that the initiation set 
point pressure was too low. 

Other output, such as the list file produed by every use 
of the ACSL was useful in the process of "zeroing in" the 
model to best predict the facility's performance. In 
particular, the DEBUG command, which provides a listing of 
most of the MMS internal variables, proved to be of great 


use in executing the models. 
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Chapter 7 


CONCLUSIONS 


The results of the previous chapter give conflicting 
impressions about the ability of the MMS to successfully 
model a small pressurized water reactor plant transient. 
However, no mention was made of the reasons why the model 
did or did not perform as desired. Because the MMS worked 
under some conditions and did not under others, a variety of 
conclusions can be drawn. The failure to operate under 


severe transient conditions will be addressed first. 


Tied Causes of MMS Failure 


4 el Major Contributions 


The MMS failed to operate past 30 real seconds when 
executing the L3-5 predictions. This failure was the most 
difficult problem encountered in this study. Determining 
the causes was the first step toward eliminating the 
problems, and proved useful in drawing conclusions about the 
MMS. 

Reference [5] indicates that if saturation conditions are 
encountered in the reactor module, the module execution 
becomes unreliable. A review of the predicted conditions at 
all module junctions in the primary system model, as well as 


those within the reactor module itself, showed that no 
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Saturation conditions were predicted before execution was 
halted. Thus an obvious cause of execution failure was 
eliminated. 

When the valve module, FRV, controlled the feed flow, a 
command to shut this valve caused the upstream pressure to 
rise to the critical point. This is the point where the 
differences between the physical states of water become 
difficult to define. Here MMS execution becomes unreliable, 
but does not terminate. Eliminating the feed flow 
regulating valve module and its associated upstream 
connection module required that the feed flow itself be 
entered as a boundary condition. This eliminated some of 
the "pure modeling" done by the MMS. It did, however, solve 
the problem of reaching the critical pressure at the inlet 
to the feed regulating valve. Unfortunately, even without 
these modules, the model still would not operate past 30 
real seconds. 

A representative of the vendor which supplied the MMS 
recommended that the selected steam generator module used in 
the L3-5 configuration be replaced by the less capable UTSG 
module, uSing the configuration of Figure 15. This step, 
too, did not allow execution past the nominal 30 seconds. 
It was noted, however, that the "low order" steam generator 
module produced very similar predictions of steam generator 


pressure as did UTSGR. 
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Finally, attempts were made to piece the experiment 
together, one separate problem at a time. First, the 
reactor was shut down. As was shown in Figure 33, the model 
would operate under these conditions. (Operation was not 
very accurate, but it continued to the desired time.) Next 
the pump coastdown was added. In this case execution 
continued to almost 40 seconds of real time. The last 
attempt added the leak, but with the pumps left in 
Operation. (ThiS was experiment L3-6 at the LOFT facility.) 
Again, execution terminated at about 30 seconds. [In the 
cases where the model would stop earlier than desired, the 
CPU time used was on the order of 25 seconds. In those 
cases where termination was at 200 seconds, CPU use was 
about 15 seconds. 

The final factor considered in examining the L3-5 
performance is the upturn in the pressurizer level in Figure 
27. This sudden change cannot be explained by any physical 
phenomena, nor is it reflected in most of the other internal 
variables. Those variables that are affected are mostly 
pressurizer module variables, which are direct inputs to the 
differential equation matrix solving subroutine of the 
FORTRAN program. Except for this sudden change in 
pressurizer states, no unusual physical properties are seen. 
The predicted trends, although not accurate, do move in the 


expected directions. 
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In all the runs which terminated earlier than expected, 
the reason noted by the CMS operating system was a FORTRAN 
error code between 245 and 281. These are the codes which 
indicate an error in the use of double (or more) precision 
variables. The FORTRAN command which overrides these types 
of errors was ignored when executing under the ACSL 
structure. 

When summed, these factors indicate that the MMS, using 
the Gear's Stiff algorithm, will continue to divide the 
differential time element until either the minimum allowed 
is reached or a FORTRAN problem with number precision is 
encountered. In none of the LOFT facility cases was 
execution ever terminated by reaching the minimum allowed 
time period. The problem is believed to be that some of the 
derivatives determined to apply to such short time intervals 
are recalculated enough times so that the storage locations 
assigned on the disk in use begin to overlap. fThis 
condition causes a FORTRAN interrupt. 

On the other hand, it should be noted that the RX3 module 
allows very rapid changes in the level of reactor power, 
implying that large derivatives over small time intervals 
are allowed. Such a rapid change is seen in Figure 29. 
This change occurs so rapidly because power level is not a 
participant in the solution of the differential equations. 
The builders of the MMS instead used the concentration of 


delayed neutron precursors and decay heat causing fission 
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products as the values that are varied by differential 
equations. The power level is computed by a single FORTRAN 
statement, using the various reactivities from the rods, 
water temperature, and boron concentration directly. Hence, 
some variables can change almost instantaneously, while 
others cannot. A determination of which variables are of 
what type is required to learn if the model has stopped due 


to excessively large derivatives. 


aomlen 2 Minor Contributions 


Problems with model parameterization were most evident in 
two recurring variables: pipe and valve flow conductances, 
and heat transfer parameters. In the case of the former, 


the conductance is typically determined by 
FC = W/p(Ap) 1/4 EG@ne yc. 


However, this equation yielded results sometimes far from 
the values eventually settled on by trial and error for use 
in the model. For example, in the hot leg piping, the 
actual pressure drop is 4 psia from the reactor outlet to 
the steam generator inlet. At a nominal density of 44 
lbm/ft and flow rate of 3.78 x 106 lbm/hr, the flow 
conductance is 4.3 x 10%. The value used in the model is 
5.0 x 10° to arrive at the same flow rate and differential 
pressure. Although this variable was the source of some 


difficulty when initially setting up the steady states, it 
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did prove useful in varying flow rates through such valves 
as BRK. 

In the case of the heat transfer coefficients, none 
proved more troublesome than those of the steam generator 
modules. After solving the series of equations provided in 
reference [5], again a trial-and-error process was needed to 
allow even the individual steam generator module to operate 
at the desired steady state. In UTSG, the final equation 


used for secondary heat transfer is 


HTC = Hbygiling/ATLMexP(Ps1/630) Eqn. 7-2 


Using as input the values of area, flow and water properties 
of references [7] and [12], this equation yields a value of 
7.8. The value needed in the L6-3 model to achieve steady 


state is 10.0. 


7.2 Satisfactory Results 


In the case of experiment L6-3, the MMS clearly shows 
that it has a capability to predict small PWR plant 
performance under some conditions. The key features proven 
useful in this effort are the automatic control functions, 
of which there are many at both the LOFT facility and full 
Sized plants, and the ability to easily change input data 
Once a model has achieved steady state. (The subscripted 
variables described in Section 7.1.1 are the exception to 


this rule.) 
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The differences in the data of Figures 29, 30, and 3l, 
while of substance, do not preclude the use of the MMS in 
performance predictions. When comparing the results of this 
study to the objectives of the MMS listed in Section 2.1, it 
is the author's conclusion that the MMS has a limited but 
reliable capability to model the thermal-hydraulic 
characteristics of a small-scale pressurized water reactor 


plant. 


73 Comparison with Other Modeling Systems 


Prediction data produced by the RELAPS/MOD1 code is 
available for experiment L3-5. This code successfully 
predicted the parameters of the LOFT facility well beyond 
the 200 second mark, So any comparison with the preformance 
of the MMS is very tenuous. In general, RELAP predicted the 
trends of the major parameters, but, similar to the MMS, at 
times the predicted and actual values were not close. 
RELAPS is a very complicated code of over 200 FORTRAN 
Subroutines, compared to the four used by the MMS. Because 
the MMS would not operate ina loss of fluid environment, 
RELAPS proved the superior under these conditions. 

RETRAN data is available for experiment L6-3. In this 
case the MMS did a better job of predicting the LOFT 
Facility performance simply because it correctly predicted 


that the reactor would shutdown automatically on low primary 
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pressure. RETRAN, although predicting that steam flow rate 
would reach a steady state of 110% rated flow, reached a 
minimum pressure of only 2117 psia, just above the scram set 
point. The RETRAN "calculated heat transfer was less than 
in the experiment, causing the calculated cooldown to be 
less severe than measured."!1 The causes of the inaccurate 
heat transfer calculated by RETRAN were not described by 
Nalezny, but it can be assumed that improper 
parameterization was a contributing factor. All codes of 
this nature suffer from this problem, including the MMS. 

No data are available on how much CPU time RELAPS5 and 
RETRAN required when performing these specific predictions. 
It is safe to say they used much more than did the MMS, 
based on studies of references [10] and [13]. 

Use of the MMS complements, rather than replaces, the 
functions of these other more sophisticated computer codes. 
In arriving at general plant design parameters the MMS 
appears to be superior because of its relatively low 
computer costs and ease of operation. The prediction of 
actual severe transient performance for in-depth safety 


anaylsis is best left to the RELAPS/RETRAN series. 


11 MNaylezny, sp. 9-19. 





Chapter 8 


RECOMMENDATIONS 


In this study, the major questions left unanswered can be 
separated into three groups: how do the volume calculations 
affect the MMS's performance?; what are the time limits to 
which a moderately severe transient can be predicted?; and 
what is the solution to the FORTRAN digit precision problem? 

Although there are thermodynamic problems with the MMS 
other than those caused by poor volume calculations, the 
effects of this parameter being determined incorrectly are 
of immediate concern. As the problem with the leak flow 
rate in experiment L3-5 showed, the total mass inventory 
does not seem to be a factor when calculating pressure 
Changes. Instead, fluid masses seem to be a module specific 
Characteristic. Further, this characteristic appears to 
apply to only the storage modules. As was noted in Chapter 
6, changing the volume of the purely resistive pipe modules 
had no effect on the rate of pressure change, while the flow 
rate out of the leak had a profound effect. Before further 
use of the MMS can be made in investigating accidents which 
involve loss of mass inventory, the method of determining 
the inventory must be corrected. The first efforts in this 
direction should be to vary the flow from a simple storage 
model, for example the stand alone pressurizer of reference 


Gu). 
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Finding the time limits of the MMS would seem to be a 
Simple problem. The difficulty, is introduced, of course, 
by also finding the limits of transient severity which the 
MMS can endure. Since at the 200 second point of experiment 
L6-3 the LOFT facility had returned to near steady state, 
the "severity limit" can be initially placed between that 
caused by a rapid increase in steam demand, and that caused 
by a small break loss of coolant. The LOFT facility's 
series of experiments is ideal for use in such a 
determination. The actual experimental results are 
available from a variety of sources in both raw and fully 
interpreted forms. Other transient predictions which the 
MMS has performed accurately on full-scale plants such as 
the Peach Bottom turbine trip can be extended to include 
those performed at the LOFT facility. Since there are 
Currently a variety of steam generator, pressurizer, 
reactor, and pump modules available, the best combination 
For each type of transient needs to be determined to utilize 
the fullest potential of the MMS. 

Investigating the problem of the FORTRAN interrupts 
requires a joint effort of nuclear and systems engineers. 
Whether or not to place a high priority on this problem is 
also a matter of question. Since there are codes now 
available which have the capabilities that the MMS showed in 
this study, solving this problem involves trading off the 


cost with the expected benefits. 
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Appendix A 


LIST OF SYMBOLS 


t - time 
é - flow entering a module (subscript) 
WwW - mass flow rate 


ij ~ Girection vectors (subscript) 


~ - energy per unit mass 

WwW - work rate 

Go ~- dimensional constant 

g - gravitational acceleration 
q - heat transfer rate 

A - heat tranfer area 

h - enthalpy 

L - length 

fy - volume integral 

> - any thermodynamic property 


A - area vector 





B 


D 


- total enthalpy 


- reverse flow 


~ effective delayed neutron fraction 


~ delayed neutron group (subscript) 


- neutrons passing between adjacent nodes 


- reactivity term 


-~ density 


- flow leaving a module (subscript) 


volume 


Smevie LOCIty 


~ heat generation per unit volume per unit time 


- shear stress 


- pressure 


- angle between flow path and a horizontal plane 


- heat transfer coefficient or internal energy 


eEomoe rh ac ure 


Seer iction factor 


20 





D - diameter 

J; - surface integral 

S - surface (subscript) 

M - mass 

Qh - partial of density with respect to enthalpy 


Setotan pressure 


On - partial of density with respect to pressure 
constant enthalpy 


Z - thermodynamic property relationship term 

12 - adjacent control volume (subscript) 

nj —- number of neutrons in node 1 

A - neutron generation time 

Cj; - delayed neutron precursors in group j and node i 
mimes mode indicator (subscript) 


r - delayed neutron group decay constant 


2 


at 


at 





122 


ieeeeeooule’s constant 
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Appendix B 


THERMAL~HYDRAULIC COMPONENTS AND NAMES 


Reactor - RXKX 
Pressurizer - PZR 
Pipes: 
hot leg = ~RXO 
cold leg - SGP 
additional hot - SSG 


leg for L6-3 


Valves: 

Soc y eee one 
pressurizer relief —- REL 
main steam control - MSS 
main stem relief - MSR 
HPIS inlet = eG 
Simulated break - BRK 
feedwater regulating SERV 

Surge junction SUE 


Reactor coolant pumps —- RCP 





Steam generator: 
experiment L3-5 


experiment L6-3 


Connections: 
coolant pump discharge 
to reactor inlet 


feedwater inlet 


HPIS junction 


Flow dividers: 
cold leg/break/spray 
steam control valve 


inlet 


Ji 


jRawr 


PRX 


REW 


ACH 


RXI 


MSL 
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Appendix C 


MMS/ACSL MODEL FILE STRUCTURE 


PROGRAM TRANS 
DYNAMIC 
LOGICAL OMASK 


CINTERVAL CINT = .1 

NSTEPS NSTP = 100000 

MAXTERVAL MAXT = 100. 

ALGORITHM IALG = 2 

CONSTANT TSTOP = 10. jeter 
OMASK = .TRUE. 

peRMt(T.GE.TSTOP) 


i] ] 
DERIVATIVE 
' START TRANSIENT CLOCK WHEN TRANS=.TRUE. ' 
LOGICAL TRANS 
CONSTANT TRANS=.FALSE. 
PROCEDURAL (ZTIME=T, TRANS) 
ZTIME=0. 
IF (TRANS) ZTIME=T 
TABLE 
Meal /0.,.2,.4,.6,.8,16,194,1.36,1.31,1.23,1.13,1./ 
Tens HVNRCP,1,8/0.,.143,.286,.429,.571,.714,.857,l.e, «we. 
-.68,-.56,-.42,-.23,--03,.2,.57,1./ 
TABLE HADRCP,1,6/-1.,-.8,--.6,--4,-.2,0.,2-54,2.03,1.82,1.61, 
1.48,1.4/ 
TABLE HVDRCP,1,8/-1.,--857,-.714,-.571,-.429,-.286,-.143,0., 


2545 DG ee Pen ae leo 1.18,1.07,.93/ 
TABLE 
eee 1. 7/0 2 4, «45 soy Op.) 225728). 33, .277.47,<71,1./ 
RaPGeE HVTRCP,1,10/0.,.11,-.22,-33,.44,.55,.66,.77,.88,1., 
BOG. Oli 89 pe eo.e370e4,.85,.89,1./ 
TABLE HARRCP,1,6/-1l. ,--8,--6 7-24, —-22,00,—le,—s6yredy cee 
.05,.13,.25/ 
TABLE HVRRCP,1,6/-1.,--8,--6,—--4,—+2,0+,~-ls 7-297,-+93, 
-.88,-.79,-.68/ 
RAPER ANREP, 1,6/0.,.2,-4,-6,-8,l-,+6,7+637+737+83,+92,1./ 
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Meme TVNRCP,1,7/0.,.2,.4,.5,-6,.8,1.,;~--48,-.36,-.26,.12,.3,; 
6 4b 
TABLE 
TADRCP,1,6/-1.,-.8,-.6,-.4,-.2,0.,2+,1.39,1.04, .8,.67,.6/ 
TABLE TVDRCP,1, ae 1. i-/9.-. oe Ce ee a eee 
a 9,1.73,1.58/1. Ee ans 5) 2G) 
TABLE TATRCP,1,4/0. eer oe 34/ 
TABLE Me i) 55" 33, 44, 55.66,.77,.89,1., ia: 
1.26,1.17,1.07, .98,.9,-78,+-67,+55,.44, .34/ 
TABLE TARRCP,1,4/-1.,-.4,-.1,0.,-l.,--91,-.52,-.48/ 
TABLE TVRRCP,1,4/-1.,-.4,-.07,0.,-1l.,--91,-.8,-.67/ 
Ter ANReP,1;/7/0.,.12,.22,.5,.7,-91,1.,0.,+85,1.09,1.02, 
Leo / 
meaeme IVNRCP,1,8/0.,.1,.2,+3,+5,+77+97l0,0.,--02,.01,.09, 
Sone SS pata les / 
TABLE IADRCP,1,10/-1.,-.9,-.7,—+6 7725770 4,72 3,72 2,7+1,0., 
-1.17,-1.23,-2.3,-2:8,-2.92,-2.68,-2., 
-1.35,-.7,0./ 
TABLE IVDRCP,1,10/-1.,-.9,-.8;—-07 77267725770 4y—03,702,02, 
-1.17,-.59,-.52,-.32,;-.19,-.1,-.03,.01,.04,.1/ 
TABLE IATRCP,1,6/0.,.2,.4,.6,-8,1. ,0.,—-3377-65,—.94, oes 
-1.2,-1.47/ 
mame IVTRCP,1,10/0.,.1,.2,-3,-4,-5,+6,+7,+85,1.,+1,-13,.15, 
.15,.12,.07,-.04,-.25,-.7,71.42/ 
TABLE IARRCP,1,6/-l.,-.8,--6;-+4,—-+.2,0.,-1.17,-.52,7+2, 
-.03,.05,.1/ 
TABLE IVRRCP,1,6/-1l.,-.8,--6,;7+4,722,0.,-1.17,-.52,;7.2, 
-.03,.05,.1/ 
Meee ENUANRCP,1,6/0.,.2,+4,-6,-8,1-,+6,.63,-73,-83,-92,1./ 
MeaeMeEOVNRCP,1,7/0.,.2,.4,-5,-6,-8,1.,-.48,-.36,—-.26,.12,.3; 


64, 1.7 
TABLE 
Week CP,1,6/-1.,-.8,-. Oe a yey. 2 eee See 4, Oy OU, 0 / 
TABLE UVDRCP,1, Bay = 1 ct) ne wo. Hee ey Soy Ol . 


pn nil oi apace 1.52,1. 38,1. )ipibe AGy/ 
TABLE UATRCP,1,4/0.,-4,+5,1.,—-68,—-27,0.,+34/ 
TABLE UVTRCP,1,10/0.,.11,.22,.33,.44;.55, .66,-77,-89,l., oe 
We 2Gp wee 7 759859, s 787007 7<55,.44,.34/ 
TABLE UARRCP,1,4/-1l.,;--4,-+1,0« ,-1ls ,--91,-252,7+48/ 
TABLE UVRRCP,1,4/-1l.,-+4,7+07,0.,—-l+,7+91,7.8,-.67/ 
MapRPeuMTROEP,1,13/0.,.05,.1,.15,.2,.3,.4,.5,+6,+7,.8,+9;,1., 
NEO eo GENO mI iGO, 73;.01,.857- «+s 
Nesp oil, .0G/ 





TABLE 


14 I 


irl ler Se, 6 Dy Ope ples, 624 Sly «337237 ~-24,0./ 
i] q 


PeMe4O( RCP’, *CPI', 'RCPO',O,0) 
q 


CONSTANT KKVRCP = 1.5E-5 
CONSTANT KNPRCP = 2. 
CONSTANT KNRRCP = 3530. 
CONSTANT KQRRCP = 5000. 
CONSTANT KPRRCP = 315. 
CONSTANT KTRRCP = 369. 
CONSTANT KVPRCP = 3.5 
CONSTANT KRRRCP = 38.31 
q 4 
CONNI('PRX','CPO','PORI') 
’ 


MeesR( ' SGP", "SGO','CPI',0,0,1) 


CONSTANT KCKSGP = .TRUE. , Keo Gres. 0G Seems. > 
KDHSGP = Q. , KLEPSOr == o> .69 yee. 
KAFSGP = .6827 , KVPSGP = 23. +7 «+. 
IHCPI = 540. 


Pier xt’, 'PORI','CLI1', 'CLI2','CLB1','*CLB2','PZSP', 'BRAK') 
q i] 


WeaLVEI(  PSY°,'PZSP’','PSP') 
CONSTANT KCPPSY=9.E66, 
KCKPSY=.TRUE., 


KCVPSY=241.4, ... 
KDHPSY=10.73 , 


KVAPSY 3 


VALVED('BRK','BRAK', 'BRST') 


CONSTANT KCKBRK -TRUE. , KCVBRK 300. , 


KKMBRK 


KVABRK = 3. 
Mri, 8/0.,25.,33-3,40.,50.,66.7775-.,100.,. 
1.,.347,.125,-0.008,-0.115,-0.103,-0.046,0./ 
Maou 8/0; 2507 33+3,40-.,50.,66.7695.,100.,... 
EOS. 8401664, .42ane0o 4, 012-07 
haere, 8/0.,252734.3,40.,50.,66.7,75.,100.,... 
leeomton 9 77.956,1.007,.778,.596,0. / 


* ¢@ 


Knee 1,3/0.,.5;eel.;1.22,1.4/ 


KRIRXX,1,3/0.,50.,100.,2258.,2258.,2258./ 
KR2RXX,1,3/0.,50.,100.,2258.,2258. ,2258./ 
KR3RXX,1,3/0.,50.,100.,2258.,2258.,2258./ 
KR4RXX,1,3/0.,50.,100.,2258.,2258.,2258./ 


Ke eeeno/ 0s, 2.,2550,42.,92.,602.0,7/6.,100.,... 
oops lO; 2 Obes tc, a0. 16,-0.13,0.12,0./ 
MEA o/0.,20.,36.0,48.,60.,74.0,86.,100.,... 
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eco lS Oo lee 2 Oeeo oO, 227 ,~-.03,0./ 
Moar KE Shxk,1,6/0.,24.,32.0,50.,67.0,100.,... 

le. 2 yO 0, Oo, Orne, 
ce  Wieneeenle, 35/0). , 50. 11002, 25521.15,25521.15,25521.15/ 
Mem KO2ZRXX,1,3/0.,50.,100.,26031.5 ,26031.5 ,26031.5 / 


ieee KLSRXX,1,370.,50.,100.,33153.9,33153.9,33153.9/ 


meme. , CLil’, CLi2’,*CLB1','*CLB2', 'HLI1','HLB1', 
'BOR',0,0,0) 


CONSTANT GBOR = 1345. 
CONSTANT KBBRXX=.000209, .001416,.001309,.002727,... 
-000935,.000314 
CONSTANT KBERXX=.000209,.001416,.001309,.002727,... 
.000935,.000314 
CONSTANT KCCRXX = .008 
CONSTANT KCMRXX = 82. 
CONSTANT KDIRXX = .43908,KD2RXX = .37078,KD3RXX=.19014 
CONSTANT KD5RXX = -2.619 , KD6RXX = -.00423 
CONSTANT KD7RXX = 3.68E-4 , KD8RXX = 8.66E-7 
CONSTANT KEPRXX = 7. ,KFCRXX = 5.4E4 
CONSTANT KFFRXX = 1.4E5 
CONSTANT KGDRXX = 3.17, 2.14 , 3.57, KGIRXX = 1.547E9 
CONSTANT KLBRXX =.0125,.0308,.114,.307,1.19,3.19 
CONSTANT KLCRXX = 5.5 , KLDRXX = 5.55E-2 , 4.3E-3 , 
6.66E-5 
CONSTANT KLERXX = .0125 , .0308 Ibe 807 7-119), 
3.19 
CONSTANT KLIRXX = 2.85E-5 , KLPRXX = 6.45 
CONSTANT KLTRXX = 11.95 , KLXRXX = 2.10E-5 
CONSTANT KMBRXX = 308. Le nO. 30 ie 
=2 = 2 OE or Or Gis ee 
Oly Or 70's 
CONSTANT KMXRXX = 1. KRPRXX = 2258. 
CONSTANT KM2RXX = -10.293 , KM3RXX = 0.0126 
CONSTANT KTSRXX = .5 
CONSTANT KT1RXX = 3.08 , KT2RXX = .07 
CONSTANT KVBRXX = 43.94,KVRRXX = 3.5, KVTRXX = 31.63 
CONSTANT KO4RXX = 0.217 , K1XRXX = 4.9117E7 
CONSTANT K10RXX = 134.01 , K1L4RXX = 2.767 
CONSTANT K2XRXX = 6.3210E-7 , K23RXX = 8.276E-4 
CONSTANT K3XRXX = -1.387E-12 , ZRORXX = 0.,0.,0. , 
Mee eel, 61.,81.,81.,81. 
‘INITIAL CONDITIONS ' 
’ § 
CONSTANT ZIPRXX = 15059 , 46106 , 2912 , 17039 , 
52156 , 3295,21043,64398,4067, 
ae BOC FOO. 7 M100. » ... 
MOG 552 , 568.4 moS.9 , ... 
WeDRxX — 5715. 95.4564, 2.E6 , 





foeeneenerigi, |,4/400.,450.,500.,920.,1.075,1.12,1. 


Ns 2. 


IHURXX = 585.1 , IPHLI1 = 2155.8 , 
TOUR KX = 2175.0 ; ZHSRXX = 540.0, ... 
ZIXRXX = 2.0E15 , 1.8E15 , 1.9E15, . 
ZIIRXX = 5.9F15, 4.9E15 , 5.5E15 
BOUNDARY CONDITIONS ' 
CONSTANT WHLB1 = 0. 
POWER = ZATRXX/1.1 
q 
PIPER('RXO','HLI1', 'HLI2',0,0,1) 
CONSTANT KCKRXO = .FALSE. , KCFRXO = 5.0E5, ... 
KDHRXO = 0. , KAFRXO = .6827 , ... 
KVPRXO = 13.56 , KLPRXO = 19.86 , ... 
IHHLI2 = 585.1 
SURINC('SUR', 'HLI2','SGI','PSG') 
CONSTANT IPHLI2=2154.8, KVTSUR=.44 
PZRB('PZR','PSG','PSP','PRF','EHTRS',3) 
CONSTANT KLHPZR = 1.25, KRCPZR = 1.417 ,... 
KLTPZR = 5. , KLRPZR = .729 , . 
KLSPZR = 6.917 , KLUPZR = 1.083 , 
KLBPZR = .708 , KAPPZR = 13.4 , ..- 
KACPZR = 6.31 , WADPZRe= 71.3 4, .<. 
KVHPZR = 3.8, KVTPZR = 39.15 ,... 
KUIPZR = 3. , KUCPZR = 110. , . 
KUDPZR = 40. , KUEPZR = 90. ., . . 
KCWPZR = .006 , KCFPZR = .56 , ... 
KCGPZR = .47 , KLLPZR = .5 , .«. 
KUTPZR = 6.27 
‘INITIAL CONDITIONS' 
CONSTANT IPPZR=2154.8 : 
ZIMPZR=838.9 ; 
ZITPZR=646.8,646.8 , 
IHPZR = 690.07,1124.00 
VALVED('REL', 'PRF', 'OUT' ) 
CONSTANT KCVREL=97.,  KVAREL=3 , KKMREL=1.0 


CONSTANT KCKREL=. TRUE. 


‘BOUNDARY CONDITIONS ' 
CONSTANT POUT=14.7 


Papeete, 2/400.,550.,.113,.113/ 


TABLE KMRITL,1,2/400.,550.,480.,480./ 


Paso tine, 2/400.75502,23.7,23.7/ 


iw yeas, 
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Memon, KTCITL,1,4/400.,450.,500.,520.,1.07E-4,1.04E-4, 


On 9E-5,9 ° Sh // 


Pome AVELTL,1,4/400.,450. ,500.,520.,9.12E-5,8.05E-5, ... 


7.17E-5,6.9E-5/ 


BmiSscR { ITL','SGI', 'SGO','FWI','STO",0,0,0) 


CONSTANT KAPITL 


KASITL = 
KCFITL = 1.385E4 , 
KDSITL = .05000 , 
RAS = «57. D> 
MRPITL = .5,.5,.5,.5% 
KPMITL = 194.2 , 
KPAITL = 
KAIITL = 
feed, 1.9,1.5,1.,1.,1.,1l. 
KLIITL = 
mime, 50,.545,.735,1 
ReGRrE = 932.2 , ss. 
KLHITL = 7.04 , 
KNDITL = 0.5 , 
KCOITL = 1.2 , 
KSLITL = l. , 
' INITIAL CONDITIONS ' 
CONSTANT ILSITL = 2.5,4.00 , 
ILDITL = 10.4 , 
PALL = 593620; 
IHSGO = 540.0 , 
MAZE = 343.5 a 
EPUITL = 808. , 
IPSGO = 2123.3, 


"HEATER AND SPRAY CONTROLLER' 


'TOTAL HEATER 
EEHTRS = 


1.626 


v 


6.21,3.63,17.63,2.5,... 


SOBRE) IE IEF 
eles J ABIL: 


KHBITL 
KRDITL 
KMSITL 


TAVIS 18) (AACA NF 23 L032) moe 


INPUT ' 


Bhatt HaeZ 


ONOFF('HT1',PPZR,EHT1) 


CONSTANT KONHT1=2235.,KVNHT1=12. 


gp Aer een, he, Ae ss 


= Se ae peter a Lee ee reed leper) Cera a ° 


-OSs5 J eee 
oDpeDpe Ie Oe 


-6,.6 y eee 
-0082 , ... 
PAN 2 er ce 


Kel =" 206. pcs el 
RVR = 1a6e “mews 
OO IE = Seow, 6 Vas 
KSTITL = .0524 
Wiebe 6.4, “2... 
PS tie 975.0, ss. 
Hae 542.5), 2 
HE Ges oS1. 6 7. .s 
IHDITL = 485.5 , ... 
PEELE = 500.,533., 030.0, 
5) B)8 Paya ° 
IWDITL = 750.00 


KOFHT1=2265.,KVFHT1=0. 


ONOFF(*HT2' , PPZR, EHT2 ) 


CONSTANT KONHT2=2230.,KVNHT2=36. 


KOFPHT2=2245.,KVFHT2=0. 


ONOFF('SVC',PPZR,CYPSY ) 


peee 


eee 








Hie 


CONSTANT KONSVC=2275.,KVNSVC=1.0,... 
KOFSVC=2250.,KVFSVC=0.0 


ACT('PSY') 
CONSTANT KATPSY=2.0, KTCPSY=1.0 


ONOFF('RVC' ,PPZR,CYREL) 
CONSTANT KONRVC=2410., KVNRVC=1.0 
CONSTANT KOFRVC=2390., KVFRVC=0.0 


ACT('REL') 
CONSTANT KATREL=2.0, KTCREL=3.0 
CONSTANT IYREL=0.0 
VALVEC('MSS','STO3','STOL') 
CONSTANT KCKMSS = .TRUE. ,KCVMSS = 1.E5 ,KVAMSS 
KVCMSS = 2000. , KXTMSS = l. 


il 
WW 
= 
e 


VALVEC('MSR','STO4', 'STO2') 
CONSTANT KCKMSR = .TRUE. ,KCVMSR = 1.E5 ,KVAMSR = 3 ,... 
KVCMSR = 2000. , KXTMSR = 1. 
ONOFF('MSC', PSTO, CYMSR) 
CONSTANT KONMSC=1000.,KVNMSC=1.0,... 
KOFMSC= 970.,KVFMSC=0.0 
§ q 
ACT('MSR') 
CONSTANT KATMSR=2.0, KTCMSR=3.0 
ONOFF('XCC',PPZR,CYXC) 
CONSTANT KONXCC=1909.,KVNXCC=1.0,... 
KOFXCC=2250.,KVFXCC=0.0 
ACT('XC') 
CONSTANT KATXC=2.0, KTCXC=3.0 


mney XCI°, 'CPO', 'HPI', 'RCPO' ) 


VALVEI('XC','HPA', 'HPI' ) 
CONSTANT KCPXC =250. , KCVXC = 250., HHPA = 70. , 
KEKAC) =. TRUES een DHAC =10.73 , KVAXC = 3 
PHPA = PCPO+20. 


Sey ( MSL’ ,‘'STO',;'STO3', ‘STO4' ) 


VALVEI('FRV','FWI1','FWI') 
CONSTANT KCPFRV=1.E4 , KCVFRV=1.E4 , ... 
KCKFRV=.TRUE., KDHFRV= 0. , KVAFRV = 3 
§ q 
CONNI('RFW','FWI2','FWI1') 
CONSTANT WFWI2 = 202500. 
aes lease 5. 2,7.2,10000.,81.,81.,0.,0 
BEE Ba lea eoi2\, 7 .2,00000.,81.,81.,0.,0 
DereEeees, 14/0, ,5.2,742,10000.,81.,81.,0.,0 
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meme BY4,1,4/0.,5.2,7.2,10000.,81.,81.,0.,0./ 

TABLE BY5,1,4/0.,5.2,7.2,10000.,81.,81.,0.,0./ 

MIPEe BREAK,1,4/0.,9.99,10.,10000.,0.,0.,1.,1./ 

aieme STEAM,1,4/0.,6.2,16.2,200.,1.,1.,0.,0./ 

TABLE FEED,1,4/0.,6.2,7-2,10000.,1.,1.,0.,0./ 

TABLE PMPSD,1,7/0.,10.8,15.8,20.8,25.8, 30.8, 39.3, ... 

3025.,3025.,2087.,1528.,1089. 1696. 40./ 

TABLE ACTFLW,1,10/0. als Die 2S 0 35 20. , 100. 
150.,10000., .. 
3.78E6,3.78E6,1.78E6,1.19E6,8.71E5,4.51E5, ... 
4.11E5,4.51E5,4.04E5,3.96E5/ 

TABLE PRZLVL,1,5/0+,6.5,11.«342,10000.,4.24,4.24,4.,0.,0./ 

Maem CLDPRS,1,10/0.,6.5,10.,14.,30.,40.,55.,87., «es 
TOW,10000., ... 
Dismiss 2103.,1777.,1450.,1170.,1088., 
MOSS. dons. , 1447/ 

iene SGPRES,1,9/0.,7.,30.,87.,98-,110.,140., 
165 P2900", ... 
808.,808.,991.,1012.,951.,972.,986.,986.,509./ 


BSGPRS = SGPRES(ZTIME) 
BPRIFL = ACTFLW(ZTIME) 
BLEVEL = PRZLVL(ZTIME) 
BCLPRS = CLDPRS(ZTIME) 


BSPD=PMPSD(ZTIME) 
BFEED=FEED(ZTIME) 
BSTEAM=STEAM( ZTIME) 
Y1=BY1(ZTIME) 
Y2=BY2(ZTIME) 
Y3=BY3(ZTIME) 
Y4=BY4(ZTIME) 
Y5S=BY5( ZTIME) 
BBREAK=BREAK ( ZTIME) 
PROCEDURAL (YFRV = BFEED) 
meRVY = BFEED 
ENDS S ‘OF PROCEDURAL (YFRV) ' 
PROCEDURAL (RFLOW=BPRIFL) 
REFLOW = BPRIFL 
END $ 'OF PROCEDURAL (RFLOW) ' 
PROCEDURAL (RLEVEL = BLEVEL) 
RLEVEL = BLEVEL 
END $ 'OF PROCEDURAL (RLEVEL) ' 
PROCEDURAL (RCLPRS = BCLPRS) 
RCLPRS = BCLPRS 
END $ ‘OF PROCEDURAL (RCLPRS) ' 
PROCEDURAL (RSGPRS = BSGPRS) 
RSGPRS = BSGPRS 
END $ 'OF PROCEDURAL (RSGPRS) ' 
PROCEDURAL (NRCP = BSPD) 
NRCP = BSPD 
END $ 'OF PROCEDURAL (NRCP) ' 
PROCEDURAL (YMSS = BSTEAM) 
YMSS = BSTEAM 
END $ ‘OF PROCEDURAL (YMSS) ' 
PROCEDURAL (YBRK = BBREAK ) 





is 


YBRK = BBREAK 
END $ 'OF PROCEDURAL (YBRK)' 
PROCEDURAL (YRX = Y1,Y¥2,Y¥3,Y4,Y5) 

YRXX(1)=Y1 

YRXX(2)=Y2 

YRXX(3)=¥3 

YRXX(4)=¥4 

YRXX(5)=¥5 
END $ 'OF PROCEDURAL (YRX)' 
INTEGER COUNT, KOUNT 
PROCEDURAL ( COUNT=) 

CONSTANT COUNT=0,KOUNT=1000 
COUNT=COUNT+1 
TERMT (COUNT.GE.KOUNT) 
END $ ' OF PROCEDURAL ' 

END $ ' OF DERIVATIVE ' 
END $ ' OF DYNAMIC ' 
END $ ' OF PROGRAM ' 
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Appendix D 


MMS/ACSL COMMAND FILE STRUCTURE 


Sele TITLE="LOFT REACTOR STEADY STATE" 
SET TCWPRN=72 , TJNITG = 1.E66,PRN=9 ,IALG=2 


em hMZRXX = -10.3 
SET ZIPRXX = 15204,50370,3247,16942,54825 eee 
3517, 2062565450, 4180 pes 
WACPI = 539.69 ’ Ker SURS- 1.85 povaler. 
KCFITL = 9.70E4 ’ ee Ce meee 2 yp ede 
IW1RCP = 1.89E6 ’ REPORT = 20/9625 | «ss 
WFWI2 = 2.205E5 ’ Poh@ Lo a2.Z209Rh9,. «6. 
RFWI2 = 51.73 ’ 1 DS eh ee ere oe 
PSTO2 = 808. ’ IPFWI1 = 809. yee 
PFWI2 = 809. ’ HFWI2 = 425. foe. 
WSTO = i220 555 ’ WFWI1l = 220500 , ... 
WFWI = 220500 ’ HFWI SZ LG) op caine 
YREL = Q. ’ WSGI aoe 7 GEO; 4 6 
YBRK = Q. ’ KVCPSY = 2000. , ... 
KVCBRK = 2000. ’ Key Orci 2000 6c euere 
KVCFRV = 2000. ’ N = 1 poe. 
PBRST = 14.7 ’ IYPSY = 0. poe: 
IYMSR = Q. ’ PPSG = (03). ros 
BHGESY = 1123.7 ’ PRED Sie = 57 213.9 4 Me 
PEeGPSY = 5.78 ’ KGPEor)—2 iz). Miles ts 
PeMIiTL = 97. ’ heat tiem ea, LlGe LO. 7. <a 
RMSiTL = 122. ’ MeSot ieee O242>O.0590) - 7. <2 
ILDITL = 10.5 ’ UES ee De 7 weg 
HAZITL = 539.17 ’ Mad beams 7 . 2 Oa 
IHSGO = 540. ’ PHA hee 0.) 22 een. 
IH1ITL = 846.752 ’ PAE 207.75 y ..- 
IHDITL = 482.47 ’ IPUITL = 808. pose. 
MEMETE = 572.9,537.81,524.27,520.25 yee 
MESGO! = 2123.3 ’ JOO EI iy 1 Sic) 0 ie ley A eae 
IPPZR = 2152. ’ GEERKR = Zl 56 7 ses 
Peer = 2153. ’ JE )S0e) IL, ie 7a s) 72 = eee 
Meoiiu =) 4.59715.,17.6,2.55 poe. 
KCVMSS = 1.E6 ’ RVCAG Saez 000.. 7, <2 
RHPA = 62.11 ' 50 ,4@ = Q. ewer 
IYXC = 0. 


PROCED NULL 
SET NDBUG=1,TSTOP =750.,CINT = .5, KOUNT = 400000 
PREPAR T,PPZR,YPSY,EEHTRS 
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OUTPUT 

T,WCPI,PPZR,WPSG, POWER, PSTO,WFWI,YXC,ZLLPZR, 'NCIOUT'=200 

START 

SET CALPLT=.T.,GRDCPL=.F.,SYMCPL=.T.,NPCCPL=50,TTLCPL=.T. 

SET XINCPL=5.,YINCPL=5. 

SET TITLE =" LOFT/ACSL SIMULATOR" 

PLOT 'XAXIS'=T,'XTAG'='(SEC)','XLO'=0.,'XHI'=750.,PPZR ,... 
'TAG'='(PSI)','LO'=2000. ,'HI'=2500.,'CHAR'='*' 

PLOT 'XAXIS'=T,'XTAG'='(SEC)','XLO'=0.,'XHI'=750.,YPSY, ... 
'TAG'='( OPEN)','LO'=0. ,'HI'=1.5,'CHAR'='*' 

PLOT 'XAXIS'=T,'XTAG'='(SEC)','XLO'=0.,'XHI'=750.,EEHTRS,... 
'TAG'='(KW)','LO'=000. ,'HI'=60. ,'CHAR'='*' 

SET NDBUG = 1 

CONTIN 

SAVE ‘IC' 

SPARESCONTINSSPARE 

STOP 

PROCED RUN 

RESTOR ‘IC’ 

REINIT 

SET ZZTICG = 0. 

SAVE ‘IC’ 

SET TRANS = .TRUE. , [ALG=2 

SET KOUNT = 6000000 , NDBUG = 1, TSTOP=20. ~ ,CINT=.05 

PREPAR ZTIME,RSGPRS,... 
PSTO,WSGO,RFLOW,RLEVEL, ZLLPZR, PCPI,RCLPRS 

OUTPUT T,PSTO,WCPI,WPSG,POWER,YBRK, ... 
WFW1,WHPI,ZLLPZR, 'NCIOUT'=40 

START 

RANGE 'ALL' 

SET CALPLT=.T.,GRDCPL=.F.,SYMCPL=.T.,NPCCPL=50,TTLCPL=.T. 

SET XINCPL=5.,YINCPL=5. 


Sele LtliTLE =" LOFT/ACSL SIMULATOR" 
PLOT 
'"XAXIS'=ZTIME, 'XTAG'='(SEC)','XLO'=0.,'XHI'=200.,ZLLPZR,... 


'TAG'='(PSI)','LO'=00. ,'HI'=6.,'CHAR'='*',RLEVEL,... 
'CHAR'='@', 'SAME', 'OVER' 
PLOT 
'XAXIS'=ZTIME, 'XTAG'='(SEC)', 'XLO'=0., 'XHI'=200.,WSGO,... 
'TAG'='(PSI)','LO'=00. ,'HI'=6.,'CHAR'='*',RFLOW,... 
'CHAR'='@', 'SAME', 'OVER' 
PLOT 
'XAXIS'=ZTIME, 'XTAG'='(SEC)','XLO'=0.,'XHI'=200.,PSTO,... 
'TAG'='(PSI)','LO'=00. ,'HI'=6.,'CHAR'='*',RSGPRS,... 
'CHAR'='@','SAME','OVER' 
PLOT 
'XAXIS'=ZTIME, 'XTAG'='(SEC)','XLO'=0.,'XHI'=200.,PCPI,... 
'TAG’="(PSI)','LO'=00. ,'HI'=6.,'CHAR'='*',RCLPRS,... 
'CHAR'='@', 'SAME', 'OVER' 
SET NDBUG=1 
SAVE 'CONT' 
CONTIN 
STOP 
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